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INTRODUCTION 


The  purpose  of  the  present  investigation  is  to  study  the  hydrodynamic 
properties  of  a hydrofoil  of  finite  span  moving  with  constant  velocity  through 
deep  water  at  a fixed  distance  beneath  the  free  surf  -cc.  The  span  of  the 
hydrofoil  is  parallel  to  the  plane  of  the  free  surface.  This  problem  has  re- 
cently attracted  a great  deal  of  attention  from  both  the  theoretical  and  the 
practical  points  of  view.  However,  the  idea  of  using  the  hydrofoil  as  an  aid 
for  locomotion  over  a water  surface  is  certainly  not  new. 

As  early  as  1898  Forlanini  in  Italy  tried  to  use  hydrofoils  for  the  pur- 
pose of  supporting  high-speed  boats.  When  such  a boat  travels  fast  enough, 
the  hydrofoil  system  installed  under  the  hull  can  lift  th'  hull  above  the  water 
surface  so  that  the  familiar  ship  bow-waves  are  replaced  by  much  weaker 
wives  due  to  the  submerged  hydrofoil  system.  Consequently,  the  resisting 
force  of  the  water  may  be  considerably  reduced,  and  the  s pe e d - powe r ratio 
miy  thereby  b^  greatly  inc  . :ased.  Another  practical  application  of  hydro- 
foils was  made  in  1911  by  Guidoni  * in  Italy  who  replaced  the  ordinary  floats 
of  a seaplane  by  hydrofoils.  The  use  of  hydrofoils  is  floats  provides  greatly 
improved  aerodynamic  characteristics  of  the  seaplane  during  flight.  However, 
the  lack  of  knowledge  of  the  fundamental  properties  and  design  data  of  hydro- 
foils has  prevented  the  use  of  hydrofoils  in  modern  applications.  The  recent 
revival  of  interest  in  the  theoretical  and  experimental  studies  of  this  problem 
is  aim^d  at  removing  these  gaps  in  basic  information. 

It  is  known  that,  for  i wing  moving  in  an  infinite  fluid  medium,  one  may 
neglect  the  influence  of  gravity  and  consider  only  the  ini  rtia  and  viscous 
effects,  the  Reynolds  condition  for  dynamical  similarity  then  holds.  Under 
this  condition,  the  nondimensional  lift  and  drag  coefficients  depend  only  on  the 
Reynolds  numbe  r and  the  geometry  of  the  body.  However,  the  situation  for  3 
hydrofoil  near  the  water  surface  is  quite  different.  The  hydrofoil  differs  from 
the  airfoil  not  only  because  of  the  possible  occurrence  of  cavitation  but  also 
through  the  strong  effect  of  the  free  water  surface.  For  the  motion  of  a 
hydrofoil  at  shallow  submergence,  one  must  consider  the  gravity  effect  be- 
cause the  wave  formations  on  the  free  water  surface  will  influence  decisively 
its  hydrodynamic  properties.  It  follows  that,  in  this  case,  the  nondimensional 
lift  and  drag  coefficients  will  be  functions  not  only  of  the  Reynolds  number  and 
its  geometry,  but  .also  of  the  Fronde’  number  and  the  cavitation  number. 
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Th^  program  of  this  report  is  as  follows:  After  a brief  survey  of  the 
ivailible  theoretical  and  experimental  information  on  the  characteristics  of 
‘lydrofoils,  the  theory  for  a hydrofoil  of  finite  span  will  be  formulated.  The 
liquid  medium  is  assumed  to  be  incompressible  and  nonviscous  and  of  infinite 
depth.  The  basic  concept  of  the  analysis  is  patterned  after  the  famous  Prandtl 
wing  theory  of  modern  aerodynamics  in  that  the  hydrofoil  of  large  aspect  ratio 
miy  be  replaced  by  a lifting  line.  The  lift  distribution  along  the  lifting  line  is 
the  same  as  the  lift  distribution,  integrated  with  respect  to  the  chord  of  thc 
hydrofoil,  along  the  span  direction.  The  induced  velocity  field  of  the  lifting 
line  is  then  calculated  by  proper  consideration  of  lift  distribution  along  the 
lifting  line,  free  water  surface  pressure  condition  and  wave  formation.  The 
"local  velocity"  so  determined  for  flow  around  each  local  section  perpendicular 
to  the  span  of  the  hydrofoil  can  be  considered  as  that  of  a two-dimensional  flow 
around  a hydrofoil  without  free  water  surface.  The  only  .additional  feature  of 
the  flow  in  this  sectional  plane  is  the  modification  of  the  geometric  angle  of 
attack,  is  defined  by  the  undisturbed  flow,  to  the  so-called  effective  angle  of 
attack  on  account  of  the  local  induced  velocity.  Thus  the  local  sectional 
characteristics  to  be  used  can  be  taken  as  those  of  a hydrofoil  section  in  two- 
dimensional  flow  without  free  water  surface  but  may  involve  civitation.  More 
precisely,  the  hydrofoil  section  at  any  location  of  the  span  has  the  same  hydro- 
dynamic  characteristics  as  if  it  were  a section  of  an  infinite  span  hydrofoil  in 
a fluid  region  of  infinite  extent  at  a geometric  angle  of  ittack  equal  to  , 

together  with  proper  modification  of  the  free  stream  Velocity.  Such  charac- 
teristics may  be  obtained  by  theory  or  by  experiment  and  should  be  taken  at 
the  same  Reynolds  number  and  cavitation  number.  With  this  separation  of 
the  three-dimensional  effects  and  the  two -dimensional  effects,  the  effects  of 
Froude  number  are  singled  out.  Thus  ::  systematic  ind  efficient  analysis  of 
the  hydrofoil  properties  can  be  made. 

The  lifting  line  theory  used  in  this  problem  is  developed  through  various 
stages  to  include  the  formulation  of  the  direct  or  indirect  problems  and  also 
the  problem  of  minimum  drag.  F inally,  detailed  calculations  of  the  lift  and 
drag  coefficients  of  a specified  hydrofoil  are  carried  out  for  the  case  of  ellip- 
tical distribution  of  circulation.  The  effects  of  the  free  water  surface  and 
the  wave  formation  are  examined  in  detail.  Needless  to  say,  the  neglect  of 
viscous  effects  will  cause  the  omission  of  the  frictional  drag  in  the  total  dr  ig 
calcul  ition.  However,  this  error  can  be  estimated  separately.  The  problem 
of  the  unsteady  motion  of  t hydrofoil  is  briefly  discussed  in  Appendix  II, 


SURVEY  OF  PREVIOUS  RESULTS 


The  invc  itiga  tion  of  the  properties  of  hydrofoils  has  been  so  far  mostly 
experimental.  Although  there  exist  some  approximate  theories,  they  do  not 
appear  to  show  the  effects  of  wave  formation  on  the  hydrodynamic  charac- 
teristics of  the  hydrofoil,  especially  when  the  submergence  is  shallow.  In 
Guidoni's  work',  it  is  claimed  that  some  of  the  advantages  of  hydrofoils  over 
conventional  ships  or  floats  are  better  lift-drag  ratio  and  less  sensitivity  to 
a rough  water  surface  so  that  their  use  may  result  in  a substantial  decrease 
in  the*  structural  Weight  of  the  hull.  On  the  other  hand,  a few  disadvantages 
are  also  mentioned;  for  example,  the  danger  to  the  hydrofoil  caused  by  drift- 
ing wood  and  seaweed.  At  high  take-off  velocities  for  the  application  to  sea- 
planes there  is  also  the  possibility  of  poor  performance  due  to  profile  cavita- 
tion. Some  early  preliminary  work  on  the  hydrofoil  problem  was  carried 
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out  by  Keldysch  ind  Liwrentjew  , Kotschin  and  Wladimirow  in  ZAHI. 

Their  studies  are,  in  a broad  sense,  only  an  extension  and  modification  of 
Lamb's  work  on  the  submerged  cylinder'*  and  Havelock's  work  on  the  sphere 
beneath  a water  surface^’  \ In  Keldysch's  work,  which  is  a special  case  of 
Kotschin's,  the  problem  of  the  two-dimensional  hydrofoil  was  solved  by  re- 
placing the  wing  by  a circular  cylinder  with  circulation.  Wladimirow  solved 
the  problem  of  hydrofoils  of  finite  span  by  replacing  the  hydrofoil  by  a horse- 
shoe vortex  of  constant  circulation  along  the  span  and  by  assuming  that  the 
free  water  surface  remains  flat  no  matter  how  close  the  hydrofoil  is  to  the 
surface.  The  calculation  then  reduces  to  a two-dimensional  problem  in  the 
Trefftz's  plane  behind  the  wing.  This  approximation  is  even  poorer  than  the 
.assumption  of  infinite  Froude  number.  In  1935,  a test  of  i single  hydrofoil 
(the  NACA  0.  0009  profile)  was  carried  out  by  Wladimirow  and  Frolow  in  the 
towing  tank  of  ZAHI*.  Experimentally  it  was  found  that  both  lift  and  drag 
coefficients  decrease  with  the  depth  of  immersion,  but  the  rate  of  decrease 
of  the  lift  coefficient  was  faster  than  that  of  drag,  especially  for  small  depths. 
The  comparison  between  theory  and  experiment  on  this  point  was  not  very 
satisfactory.  About  the  same  time,  the  problem  of  cavitation  on  a submergei 
obstacle  of  hydrofoil  section  such  as  a ship  propeller  was  approached  by 
Ackeret  , Walchner  , Lerbs  , Martyrer  , Gutsche  , and  Orr.ith 

Around  1937,  further  contributions  to  the  hydrofoil  problem  were  made  by 

14  15  1 6 

Wcinig  , Tietjcns  , and  von  Schertel  . Weinig  gave  a preliminary  yet 
exhaustive  discussion  on  hydrofoil  and  planing  problems  neglecting  completely 
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the  effects  of  gravity  and  the  free  surface  elevation;  the  wave  resistance  for 
both  cases  was  estimated  approximately  afterwards.  During  the  last  World 
War,  a series  of  tests  on  a number  of  hydrofoils  of  different  foil  sections  was 
carried  out  at  NACA  by  Land,  Benson  and  Ward  ’ * 7 in  the  NACA  towing 

tank.  Several  interesting  results  were  found.  At  depths  greater  than  4 
chords,  the  influence  of  the  free  surface  is  negligibly  small.  In  the  range  of 
depths  between  4 chords  and  approximately  1 / 2 chord,  lift  and  drag  co- 
efficients decrease,  and  the  cavitation  speed  increases,  with  decrease  in  depth 
until  the  hydrofoil  approaches  and  breaks  through  the  surface  at  which  point  a 
sudden  decrease  in  lift  occurs.  However,  the  corresponding  values  of  the 
lift-drag  ratio  increase  to  maximum  values  when  the  hydrofoils  are  near  the 
surface  then  decrease  rapidly  with  further  decrease  in  depth  to  values  for 
the  planing  surfaces.  Obviously,  these  observations  indicate  that  there  exists 
an  optimum  value  of  the  depth  for  which  the  hydrofoil  operates  most  efficiently. 
The  effect  of  an  increase  of  speed  is  to  reduce  the  minimum  drag  coefficient 
and  to  make  it  occur  at  a corresponding  higher  value  of  the  lift  coefficient 
until  cavitation  takes  place.  Any  lower  surface  cavitation  reduces  the  lift; 
while  complete  upper  surface  cavitation  prevents  further  increase  of  lift. 

Under  such  condition,  however,  a maximum  value  of  lift  as  high  as  1 ton  per 
square  foot  was  recorded.  A profile  with  sharp  leading  edge  seems  desirable 
for  reducing  cavitation.  The  hydrofoil  section  NACA  16-509  was  observed  to 
have  favorable  hydrodynamic  properties.  The  stability  problem  of  a system 
of  hydrofoils  was  studied  by  Imlay^.  An  exhaustive  exper in.ental  investiga- 
tion of  a single*  hydrofoil  with  large  dihedral  was  made  by  Sottorf^  in  the 
Gottingen  towing  tank  in  which  eighteen  different  profiles  were  studied;  his  re- 
sults have  been  made  available  only  recently.  He  found  that  a thin  profile  of 
almost  circular  segment  form  with  pointed  nose  and  an  increased  convex  camber 
in  the  nose  area  is  favorable  for  laminar  flow,  and  consequently  also  reduces 
cavitation,  gives  higher  lift  and  lift-drag  ratio,  and  minimum  spray  formation 
at  the  junction  of  the  tip  of  the  hydrofoil  and  the  free  water  surface.  The  main 
feature  of  the  depth  effects  on  hydrodynamic  properties  are  in  good  agreement 
with  NACA  results.  It  was  also  found  that  a partial  cavitation  on  the  upper 
surface  has  a favorable  effect  provided  that  it  covers  less  than  a half  chord 
because  a thin  cavitation  bubble  layer  reinforces  the  flow  curvature  and  thereby 
increases  the  lift.  In  addition,  the  bubble  layer  lessens  the  surface  friction  by 
acting  as  a cushion  between  the  solid  surface  and  the  high  velocity  water  stream. 
However,  this  favorable  range  of  operation  is  not  very  stable. 
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general  formulation  of  the  problem 

Consider  a hydrofoil  of  span  2b,  with  arbitrary  profile  and  plan  form, 
moving  with  constant  forward  velocity  U through  deep  water  at  fixed  depth 
of  immersion  h which  is  measured  from  its  trailing  edge  to  the  mean  free 
water  surface.  If  we  choose  a coordinate  system  fixed  with  respect  to  the 
hydrofoil,  then  in  this  system  the  flow  picture  would  appear  to  be  stationary 
with  a uniform  free  stream  velocity  U approaching  the  hydrofoil.  Let  the 
x-axis  be  parallel  to  the  direction  of  the  free  stream,  the  z-axis  point  up- 
ward, with  z = 0 representing  the  undisturbed  free  surface  and  the  trailing 
edge  of  the  hydrofoil  lying  between  (0,  - b,  - h)  and  (0,  b,  - h)  as  shown 
in  Fig.  1 . 


Fig.  1 - Coordinate  system  for  th>.  hydrofoil  motion. 

The  elevation  of  the  disturbed  free  surface  caused  by  the  hydrofoil  is  denoted 
by  £(x,y)  measured  from  z = 0.  The  liquid  medium  is  assumed  to  bv  in- 
compressible and  nonviscous  so  that  the  condition  of  lr rotationaiity  and  con- 
tinuity implies  that  the  perturbation  velocity  potential  satisfies  the 

Laplace  equation 

r7Z  if  = 0 . 


0) 
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From  its  solution  the  total  velocity  field  of  the  liquid  flow  can  be  obtained  as 


O - U i + q , q = (q^  q2>  q3)  = grad  <p  (2) 

The  pressure  field  may  be  determined  from  the-  Bernoulli  integral 

'£+-jC>Q  + gz  = const  (3) 

in  which  the  effect  of  gravity  has  been  taken  into  account. 

Now  let  us  assume  that  the  resulting  motion  is  such  that  almost  every- 
where I q ) is  much  less  than  U so  that  Eq.  (3)  may  be  linearized  to  yield: 


+ Ue  + bz  - const, 

p x ° 


(4) 


The  corresponding  line  uized  boundary  conditions  on  the  disturbed,  free 
surface  now  become 


and 


(5) 


(&) 


Eq.  (5)  represents  the  line  irizcd  kinematic  condition  on  the  free  surface, 
while  Eq.  (b)  is  identical  with  the  dynamic  condition  p = pQ  - constant 
on  the  disturbed  surface.  It  may  be  noted  that  gravity  affe  cts  the  kinematics 
through  the  boundary  conditions  (5)  and  (6)  and  the  r,  after  influences  the 
dynamics  through  Eq.  (*1),  even  though  it  docs  not  appear  explicitly  in  the 
kinematical  Eq.  (1).  The  boundary  condition  on  the  hydrofoil  surface  is 

rgr  - - U cos(n,  x)  (7) 

on 

where  ^ is  the  normal  component  of  the  velocity  at  the  surface.  It  may 
also  be  mentioned  here  th  it  the  real  physical  situation  requir,  s that  the  dis- 
turbance upstream  should  diminish  at  a rate  which  is  made  to  agree  with 
observation. 

In  order  to  simplify  the  problem,  it  will  be  assumed  that  the  aspect 
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of  the  hydrofoil  with  plan  area  3 and  span  2b  is  so  large  that  the  whole 
hydrofoil  may  be  replaced  by  a lifting  line.  Then  the  subsequent  analysis 
can  be  made  similar  to  the  Prandtl  wing  theory  of  modern  aerodynamics 
by  introducing  appropriate  modifications  of  Prandtl's  original  concept.  The 
fundamental  concept  is  that  the  circulation  distribution  along  the  lifting  line 
which  replaces  the  hydrofoil  is  the  same  as  the  distribution  of  bound  vorticity 
integrated  with  respect  to  the  chord  of  the  hydrofoil  along  the  span  direction. 
We  then  calculate  the  induced  velocity  field  of  this  circulation  distribution  on 
the  lifting  line.  The  characteristic  length  of  this  induced  velocity  field  is 
the  span  2b  or  the  immersion  depth  h,  whichever  is  the  smaller.  There- 
fore if  h is  much  larger  than  the  representative  chord  of  tne  hydrofoil,  then 
the  characteristic  length  of  the  induced  velocity  field  is  very  much  larger  that 
the  characteristic  length  of  the  local  velocity  field  which  is  associated  with 
the  bound  vorticity  and  the  circulation.  Then  as  far  as  the  lift  production  is 
concerned,  each  section  of  the  hydrofoil  is  effectively  surrounded  by  a 
stream  of  infinite  extent  moving  with  an  effective  velocity  which  is  not  the 
free  stream  velocity  U but  the  sum  of  U and  the  induced  velocity.  V.’e 
then  replace  the  boundary  condition  of  Eq.  (7)  by  the  statement  that  the  cir- 
culation of  the  lifting  line  at  a certain  spanwise  station  is  the  same  as  that 
corresponding  to  a two-dimensional  flow  without  free  surface  around  the 
same  local  section  with  an  effective  free  stream  velocity  equal  to  the  sum  of 
U and  the  induced  velocity.  Although  our  argument  for  the  lifting  line  con- 
cept is  based  upon  large  aspect  ratio  SR.  Experience  with  airfoils  indicates 
that  the  theory  can  be  expected  to  be  sufficiently  accurate  even  for  .''R  as 
low  as  four. 

We  shall  assume  that  along  the  lifting  line  the  distribution  of  the  cir- 
culation strength  P (y)  is  known  for  given  angle  of  attack  and  immersion 
depth.  It  should  be  pointed  out  here  that  for  shallow  submergence  the  effects 
of  the  free  surface  modify  considerably  the  approaching  flow  velocities  and 
thereby  make  the  pressure  and  lift  distribution  different  from  its  correspond- 
ing aerodynamic  value.  Therefore,  in  our  hydrofoil  problem  the  lift  distri- 
bution is  no  longer  simply  proportional  to  P(y)  as  it  is  in  aerodynamic  wing 
theory  but  is  rather  a complicate  function  of  n(y).  When  P (y)  is  given, 
the  lift  distribution  can  be  determined  as  will  be  shown  later.  However,  the 
calculation  of  P(y)  for  known  lift  distribution  is  more  difficult.  Strictly 
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spe iking,  T (y)  depends  not  only  on  the  flow  velocity,  angle  of  attack  and  the 
geometry  cf  the  hydrofoil,  but  also  on  the  depth  of  submergence  h.  When  h 
• nds  to  it  ;i  itv  so  1.  i.  *.h:a  hydrofoil  is  in  a flow  of  infinite  extent,  the  free 
••urine*  n.-vl  v.  e ..  r:y  should  vanish  and  P (y)  then  tends  to  its  conventional 
•'»>  roeyna.i'ic  value  ■ ).  5. nee  the  circulation  for  equal  airfoil  and  equal 

effective  angle  o ett.u  k . a proportional  to  the  approaching  relative  flow  velo- 
city, we  have 

:..(0.  y.  -h) 

r (y)  --  1 - rjy).  (9) 


If  I (y',  instead  of  P(y)  at  d pth  h,  were  given,  Eq.  (9)  would  lead  to  an 
integral  equation  for  T (y)  as  v/il'  be  shewn  later*  A similar  situation  results 
if  one  intends  to  find  P(y)  at  h ior  given  geometry  of  the  hydrofoil.  There- 
fore we  shall  first  consider  the-  simple  case  in  which  P (y)  is  assumed  to  be 
given. 


In  this  case,  the  contribution  of  the  lifling  line  to  the  perturbation  poten- 
tial can  be  shown  to  have  the  following  integral  representation  (cf.  Ref.  22, 
also  cf.  Appendix  !): 


= — ?r.t^Ltl  f “Vfoldr, 

2"  j-oo 


00 


sin  X x 


-oo 


d X 


cos  pi(y  -q  ) 


* r”p(.Jdn 

vJ-C  o 


roo 


cos  tr(y  - n)e 


- H 


■ + h| 


d p . 


C-3) 


It  may  be  noted  that  the  above  expression  satisfies  Eq.  (1)  -nd  has  a singularity 
along  the  lifting  line,  a jump  of  ^ (y)  across  the  vortex  sh>.*t  on  the  down, 
stream  sid*  of  the  lifting  line  and  is  regular  elsewhere  inside  the  water.  Need- 
less to  say,  vj  has  no  meaning  above  th>  free  surface  (z>£).  It  may  also  be 
remarked  here  that  sorru  other  useful  representations  for  can  be-  written 

as  follows  (cf.  Appendix  I)  : 


sig^  ( z + h) 
4'tt 


oo 


r(ri)d 


-oo 


')[  dj  c'klz  + hl  Jo(k  'J(x-fjr+(y-n)2  )kdk(lOa) 
J o J o 


-V- 


or. 


*.  ■ ^ 


r cd 


-oo 


ial 


7 — “ 2"  (*  + . ,■■  =====  = ■■==)  dn,  (10b) 

' ‘ ' ' L ' ' + (y  -r)  )2  + (z  + h)  ^ / 


( y ~ r) ) +(z  + h)< 


whv.iv  denotes  the  Bessel  function  of  the  first  kind. 

It  is  convenient  to  decompose  ip  into  three  parts 


0 = ^ ♦ *2  ♦ <3 


(11) 


In  such  a way  that  the  boundary  condition  of  Eq.  (6)  is  decomposed  is  follows 


d<f1 

3 x + 3x 


= 0 at  z = 0 , 


and 


U 


3 

3 x 


'z  - 0 


(12) 

(13) 


Physically  this  decomposition  means  that  V3  represents  the  potential  due  to 
wave  formation  so  that  «p  is  the  only  part  of  <p  which  involves  the  effect 
of  gravity,  and  it  vanishes  when  the  gravitational  effect  is  neglected.  In  this 
case  the  bound  iry  condition  Eq.  (6)  is  also  relaxed.  It  does  not  follow, 
however,  that  the  disturbance  C should  also  vanish  because  X.  is  also  in- 
fluenced by  ^ ^ ind  through  Eq.  (5).  The  remaining  part  of  (p  , 
namely,  ^ j + <P^  , equivalent  to  the  potential  of  1 biplane  system  in 

infinite  flow  with  the  upper  lifting  line  of  an  equal  circulation  F (y),  rotating 
in  the  same  sense,  distributed  along  the  image  points  (z  = + h)  of  the  real 
wing.  In  other  words,  ip^  represents  the  correction  to  <p  due  to  the  effect 
of  the  mean  free  surface  z = 0.  According  to  this  reflection,  0,  may  be 

4* 

written  down  directly  from  Eq.  (10)  by  replacing  h by  ( - li).  In  the  r ange 
of  present  interest,  we  have 
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r>oo 


cos  p (y -r) ) e^Z  ^d^  for  z »*"  0 , 


(14) 
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which  may  also  be  converted  into  expressions  similar  to  Eqs.  (10a)  and  v.  10b). 
It  may  be  noted  that  this  expression  is  regular  everywhere  inside  the  water. 


Now  it  seems  plausible  to  assume  that  v ^ and  K may  be  represented 


by  the  following  expressions: 

poo  r co 
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P (X,  p)  cos  X x d X 
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cos  p(y  -13)  e 


- g x + p 


dp 


where  a(X,  p)  and  p(X,  p)  arc  two  functions  of  p and  X to  be  determined 
by  using  boundary  conditions  (5)  and  (13).  Substituting  these  values  into 
Eqs.  (5)  and  (13),  wc  obtain: 


a(X,  p) 


V, 

1 


A 


and 

P(k.  p)  - * ~ a(k.  h)  • 

These  integral  representations  of  ^ and  4 can  be  converted  into  forms  con- 
venient for  calculation  by  introducing  the  new  variables  k and  0 such  that: 


X = k cos  0 , 


= k sin  0 
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(15) 
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rrU 
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f-l  k cos(k  x co a O)cos(k(y-Tl)sip  9)  ^ ^ 

k -it  sec  ^9 

(16) 


It  can  be  seen  that  this  value  of  4 given  by  Eq.  (16)  is  even  in  < so 
thit  the  elevation  upstream  is  merely  the  reflection  of  that  downstream. 
Therefore  it  seems  necessary  to  investigate  the  behavior  of  ' for  large  nega 
live  v ilues  of  >.  . n oruer  to  make  a comparison  with  observation.  The  in- 

tegral representing  £ as  given  by  £q.  (16)  is  indeterminate,  but  it  is 
f lrnihar  in  pi  me  wave  problems  that  its  principal  value  can  be  obtained  by 
considering  .<  to  be  complex  and  evaluating  it  by  contour  integration.  The 
path  on  the  positive  real  axis  should  be  indented  iround  the  point 
k - V.  sec~Q  uhich  is  a simple  pole  of  the  integrand.  . fter  deforming  the 
contour  to  the  imaginary  axis,  we  obtain 


i-nncipal  value  of 


.oo 


-hK  cos(k  x cost*)  cos(k(y-g)  sin  8)  ^ 

k - r.  sec ^0 


- ( si  gn  x)  it  y se  c “w  e SlC  ^sin(K.xsec9)cos(.(y-r))sec^Qsin9) 


oo 


■ k i x i co 
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s0  i/w  \ n\  ' sccu0  cosfk  h)  + k sin(kh)  ... 
cos  h(k(y-". ) sin  9)  ■*  ■ * *—  -L  k d k 

k 4 < sec  9 


(17) 


The  integration  of  the  last  integral  in  Eq.  (IT)  with  respect  to  9 (cf.  Eq. 

(16))  should  only  c-  tend  over  the  range  for  which  x cosO-f  (y-n)sinO  > 0 
for  an  assigned  t>omt  (x,y),  (cf.  App.  1!  for  the  detail).  It  is  evident  from 
Eq.  (IT)  that  "'ltiongh  '.he  integral  vanish.es  as  x • eo,  th,_  first  ttrmrcm.v.ns 
even  for  x - - o . Then  there  aie  surface  waves  even  far  ahead  of  the  hydro 

foil.  Tins  is  m contradiction  with  experience.  The  paradox  can  be  resolved 
by  observing  tii  it  there  ~ re  solutions  of  the  basic  flow  equation,  the  free  wave 
solutions,  whim  s itisfy  the  surface  pressure  condition  of  Eq.  (6),  V'e  can 
then  'dd  these  solutions  to  our  velocity  potential  without  rendering  the  solition 
so  far  obtained  invalid  for  our  problem.  The  appropriate  free  wave  solution 
to  be  added  is  determined  by  the  condition  that  the  free  surface  elevation  i 
must  vanish  at  large  distances  ahead  of  the  hydrofoil.  Therefore,  we  have  tc 
superimpose  on  £ another  system  of  waves 
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- co  r,ff/2 
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-CD  v. 


sec^O  o S"'C  ®sin(«*’.x  sec  9)  co»(/  (y -q)  sec^O  sin  0)  d 0 

(18) 


which  is  odd  in  x.  In  order  to  remain  consistent  with  all  the  boundary  con- 
ditions, we  must  add  to  the  velocity  potential  another  term,  say,  <p^,  given  by 


r>  oo 
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.•c^0e'^Z  ^Sec  ®cos(Kx  sec  8}  cos(,((y-r))  sec^O  sin  0)d  0, 

( for  z ' 0)  (19) 


which  obviously  satisfies  Eq.  (1).  This  part  of  the  potential  may  be  thought 
of  as  a correction  term.  The  same  result  can  be  obtained  by  introducing  a 
fictitious  viscous  term  to  the  Bernoulli  equation  (cf.  Appendix  II,  Lq.(  18)  ),  a 
device  first  discovorec  by  Lord  Rayleigh.  The  physical  significance  of  this 
viscous  term  is,  however,  not  easy  to  understand.  Our  argument,  even  if 
somewhat  lengthy,  has  its  merit. 

I 

From  the  resulting  expression  of  the  surface  elevation,  4 + 4 , it  can 
be  seen  that  the  term  with  the  factor  c*ktxicos0  jn  £qt  (17)  diminishes  ex- 
ponentially with  increasing  ixi  both  m upstream  and  downstream  direction. 
Hence  for  large  positive  values  of  x,  the  surface  elevation  can  be  approximated 
by 


xOo 


2 \ 

7TT 


r/z  20 

(■Hjdq  sec  Oe  sec  $in(<<  x sec  0)  cos(,\(y-q  )sec  Qsin0)d0  , 


j-oo 


( for  x > 0, large)  (20) 


This  relation  snows  that  4 is  analyzed  into  components  of  simple  waves, 
where  each  train  of  waves  (corresponding  to  one  value  of  0)  propagates  on 
the  downstream  side  of  the  wave  front  x cos  0 + | y-q  | sinC  = constant 
(cf.  Appendix  II)  with  wave  length 
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Z TT 


- sec  0 


2 tt  U'~ 
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(21) 


This  configuration  of  surface  waves  resembles  that  of  typical  ship  waves 
(cf.  Ref.  5,  pp.  433-4  17).  To  study  the  behavior  of  4 directly  behind  the 
hydrofoil  as  x -*  r oo,  y finite,  we  can  approximate  the  integral  in  Eq.  (20) 
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by  applying  the  method  of  stationary  phase  (e.  g.  cf.  Ref.  23,  p.  505).  Only 
the  result  will  be  given  here: 

-,oo 

P ("H  ) d n + C ( ~ ) , as  x -*  + oo,  y finite  . ( 22) 

-oo 

Therefore  ^ tends  to  zero  like  0(x  ^/^)  (ar  downstream  ; and  at  fixed  x 
the  effect  of  free  stream  velocity  is  such  that  £ is  proportional  to  U”*, 
noting  that  in  Eq.  (22)  T (y)  is  proportional  to  U.  Eq.  (22)  shows  that  as 
the  distance  from  the  hydrofoil  is  increased  only  the  total  circulation  on  the 
lifting  line  matters,  the  detail  of  the  distribution  of  circulation  is  mconsequei 
tial.  This  is,  of  course,  what  we  would  expect  from  general  principles. 

It  may  be  also  be  seen  from  Eqs.  (16)  and  (18)  that  there  is  a local  dis- 
turbance immediately  above  the  obstacle.  As  we  do  not  intend  to  examine 
the  surface  wave  pattern  in  detail,  the  above  discussion  suffices  to  describe 
the  general  behavior  of  the  surface  elevation. 

Thus  far  wc  have  found  the  complete  perturbation  potential,  namely, 

: Sj  ♦ c2  + V'3  ♦ <r4  (2  3) 

where  i?| , ^ and  ,*4  arc  given  by  Eqs.  (10),  (14),  (15)  and  (19) 

respectively. 
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CALCULATION  OF  THE  DRAG  FROM  INDUCED  VELOCITIES 


The  total  drag,  neglecting  the  part  caused  by  frictional  effects,  ex- 
perienced by  the  hydrofoil  may  be  calculated  by  using  the  Kutta- Joukowsky 
law; 

- oo 

D = - p f (y)  -frJdy  (24) 

oZ 

J - oo 


where  - o (0,  v,  -h)  represents  the  total  induced  downwash  at  the  trailing 


edge  of  the  hydrofoil.  If  we  split  D into  components,  each  of  which  corre 
sponds  to  the  respective  component  of  0 , we  have 


D = D j + + D ^ , 
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whe  re 


°n  = - P 


oo 


i'(y)  <rn  (0,  y,  -h)dy  , (n  = 1,  2,  3,  4) 


(25) 


-oo 


The  calculation  of  D^,  the  induced  drag  due  to  the  trailing  vortices  of  the 
hydrofoil  itself,  is  familiar  in  three-dimensional  aerodynamic  wing  theory 
and  ' s given  by  (cf.  Ref.  22): 
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(27) 


-oo 


f((w)  and  g(p)  are  thus  merely  the  Four  ic  r - coefficients  of  the  circulation  dis- 
tribution P(y)  . 

3^n 

For  the  rest  of  D,  we  first  calculate  — — from  Eqs,  (14),  (15)  and 
(19)  and  obtain: 
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jubstitution  of  thes^.  values  into  Eq.  (25)  leads  to: 
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rDi  . > h ti  1 Z > 1 

°2  + °3  = * ~i\  c [<(p)  * 8(h)  ! 


(28) 


:nd 
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D4  = TT  P r 
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-2hr<seC  0 i sec^O  sin  ©)^  + g(A  sec^8  sin  8)^  J sec^CdO 

(29) 


where  f and  g are  defined  by  L.q.  (27).  The  combination  of  and 

represents  the  total  contribution  of  the  mean  free  surface  effect,  which 
favorably  decreases  the  tot  1 dr"-g  especially  when  h is  small.  As  h tends 
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to  zero,  cancels  D^.  The  part  represents  the  wave  drag  due 

to  the  downstream  wave  formation  which  results  from  the  gravity  effect.  It 
may  be  remarked  here  that  the  wave  drag  can  also  be  obtained  by  the  method 
of  traveling  surface  pressure  (cf.  Appendix  III). 

From  the  above  result  we  jee  that  the  drag  is  expressed  in  terms  of  th^ 
Fourier  coefficients,  (f(p)  and  g(^),  of  the  circulation  distribution  P (y). 
Therefore  for  a given  ( y ) , the  total  drag  is  completely  determined.  In 

most  cases,  the  wing  form  is  symmetric  with  respect  to  the  central  plane 
y = 0,  that  is, 


r'(y)  = H-y), 


(30) 


then  it  follows  from  Eq.  (27)  that  g vanishes  identically  under  this  condition. 


By  substituting  the  relations  of  Eq.  (27)  into  Eqs.  (26),  (28)  and  (29), 
the  drags  can  be  expressed  directly  in  terms  of  the  circulation  distribution. 
Thus 
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Eq.  (26a)  is  the  familiar  drag  formula  for  the  Prandtl  lifting  line  theory. 

Eq.  (29b)  shows  that  G(q-t)')  is  symmetrical  with  respect  to  y)  and 

It  is  seen  that  the  drag  components  D^,  and  arc  independent 

of  the  gravitational  effect.  They  are  thus  properly  called  the  components  of 
the  induced  drag.  All  the  drag  caused  by  the  presence  of  surface  waves  is 
represented  by  'O^.  is  thus  the  wave  drag.  Our  calculation  does  not 

include,  however,  the  skin  friction  drag  produced  by  the  viscous  shearing 
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stress  on  the  surface  of  the  hydrofoil.  This  skin  friction  drag  should  be 
added  to  the  drag  calculated  in  this  section  to  obtain  the  total  drag  of  the  hydro- 
foil. 


CALCULATION  OF  THE  TOTAL  LIFT 


It  is  well  known  that  the  lift  exerted  on  a wing  flying  in  an  infinite  fluid 
is  given  by 


L = pU 
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P 00 
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(31) 


In  the  hydrofoil  problem,  however,  the  approaching  stream  velocity  is  influ- 
enced by  effects  of  the  free  surface  and  wave  formations.  Hence,  the  total 
lift  on  the  hydrofoil  will  differ  from  its  aerodynamic  value  L^  by  an  amount 
AL  such  that  the  total  lift 
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From  the  expression  of  0 we  find  that  Cj  and  have  no  contribution  1 1 the 

value  of  at  the  hydrofoil.  The  final  result  is 
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Subs ti tmtion  of  this  value  into  Eq.  (33)  yields: 
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where  f and  g are  again  given  by  Eq.  (27)  and  g vanishes  for  symmetri- 
cal wings.  Equation  (35)  represents  the  contribution  of  free  surface  effects, 
which  tend  to  lessen  the  total  lift. 

By  substituting  the  relations  of  Eq.  (34)  into  Eq.  (33),  a direct  ex- 
pression of  AL  in  terms  of  the  circulation  distribution  P (y)  is  obtained: 


p oo  n CD 


AL  = * ■£ 


r(n)r(ri')  Ftn-r,')  dn.dV 


(35a) 


■on  O 


-CD 


where 


p oo 


F(r)-r)')  = "7* 


2 h 


p w/2 


*dp 


cos  [h-Vip  sino]  do  (35b) 


p - i\.  sec  0 


Therefore  the  magnitude  of  AL  is  in  general  equal  to  the  magnitude  of  drag, 
i.o.,  of  one  order  smaller  thin  the  lift  L. 


In  connection  with  the  present  discussion,  we  examine  the  effect  of  im- 
mersion depth  on  the  value  of  P(y).  Substituting  Eq.  (34)  into  Eq.  (9),  we 
obtain  the  following  relation  for  a symmetrical  hydrofoil, 
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(36) 


where  P (y)  is  the  value  of  ' (y)  f-r  the  same  wing  and  same  effective  angle 
of  -attack  as  h — » oo.  If  (y)  at  depth  h is  given,  then  the  aerodynamic 
value  of  the  circulation  of  the  same  wing,  (“^(y),  can  calculated  from  Eq. 
(36).  On  the  other  hand,  if  only  P (y)  is  known,  then  Eq.  (36)  provides  an 

OO 

integral  equation  for  P(y)  since  f also  depends  on  P(y).  However,  for 
moderate  values  of  h,  (y)  and  i-*  (y)  are  approximately  equal  as  will  be 
shown  later. 


GEOMETRY  OF  THE  HYDROFOIL;  EFFECTIVE  ANGLE  OF  ATTACK 

Take  a strip  of  the  hydrofoil  of  width  (dy)  in  the  spanwise  direction  and 
with  chord  c(y)  located  at  y.  According  to  the  basic  concepts  of  the  lifting 
line  theory  explained  in  the  section  "General  Formulation  of  the  Problem", 
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the  circulation  r(y)  around  the  lifting  line  at  y is  the  same  as  the  circulation 
around  the  corresponding  hydrofoil  section.  It  then  follows  that  the  lift  on  this 
lifting  line  segment  is  the  same  as  the  lift  on  this  hydrofoil  section,  integrated 
with  respect  to  the  chord  of  the  hydrofoil.  The  local  flo  r around  this  hydrofoil 
section  in  a plane  y = constant  may  be  then  considered  as  a two-dimensional 
flow  around  the  same  hydrofoil  section,  with  free  stream  velocity  (U  + u(y),  w(y)) 
and  without  free  water  surface.  Therefore  the  local  lift  coefficient  (y)  and 
drag  coefficient  C^(y)  per  unit  chord  length  can  be  defined  by 


p [U  + u(y)  )2  c(y)  c^(y)  dy  = p [u  ♦ u(  y)J  T(y)  d y 

y pi  U + u(y)j2  c(y)Cd(y)  dy  = p w(y)  P(y)dy 

where  u(y)  and  w(y)  denote  the  induced  velocities  and  respec- 

tively at  the  hydrofoil.  The  above  relations  then  give 


and 


"(y)  = 4-iu*u(y)j  c(v)  Cj(y) 


(37) 

(38) 


where  ( (y ) is  defined  by  the  above  equation  as  the  downwash  angle.  Eq.  (38) 
indicates  an  important  feature  of  the  local  flow  that  the  actual  absolute  angle  of 
attack  a , measured  fr^m  the  free  stream  to  the  icro-mt  direction,  is  modi- 
fied  to  an  effective  angle  of  attack  a>(  measured  from  the  approaching  stream 
to  the  zero-lift  direction  on  account  of  the  local  perturbed  velocity  as  shown  in 
Fig.  Z.  The  relation  between  and  >s  given  by 

a,  = a - « (y)  (39) 

c a 


D 


hi 
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In  this  way  we  can  separate  the  three-dimensional  effects  and  two-dimensional 

effects  by  taking  the  local  sectional  characteristics  of  the  hydrofoil  at  depth 

h and  at  absolute  incidence  a to  be  the  same  characteristics  of  the  section 

a 

of  a two-dimensional  hydrofoil  submerged  at  infinite  depth  and  sustained  at 

an  absolute  incidence  equal  to  o = a (y).  Referring  to  these  two  in- 

c ^ 

cidenccs,  C^(y)  is  approximately  proportional  to  qr  or  for  most  of 
the  usual  profiles  at  small  incidences  (Ref.  17  and  21).  Moreover,  when 
vanishes  so  do  qv  and  € (y).  Hence  we  may  put 

C.  (y)  = a a = a a . (40) 

I ' ’ ' a a e e 


In  general,  all  a^,  o^,  a and  a may  be  functions  of  y.  In  particular, 

a is  constant  for  wings  with  no  geometric  twist,  and  in  addition  with  the 
same  profile  along  the  span,  a is  constant.  Combining  Eqs.  (37) -(40), 
we  have 


r (y)  = 


U c(y ) ac 

2 


(41) 


For  given  '"(y),  u(y)  and  w(y)  are  determined,  then  this  relation  gives  the 
value  of  the  chord  length  c(y)  except  for  a proportionality  constant.  In 
particul  ar. 


U 


r (0) 


(41a) 


which  gives  the  relation  between  ■~'(0)  and  cq,  the  chord  at  the  central 
s«  ction. 


FORMULATION  OF  PROBLEM  WITH  SPECIFIED  GEOMETRY 


The  problem  is  for  given  geometry  of  the  hydrofoil,  that  is,  given  b, 
c(y),  q (y),  a and  h,  to  find  P(y),  C.  and  Cn  . The  relation  between 

1 C L U 

the  g iven  quantities  and  the  unknown  f(y)  was  approximated  in  the  section, 

"Geometry  of  the  Hydrofoil;  Effective  Angle  of  Attack"  , to  be 


T 

7 


ZhL 


U c(y ) a 


v ^ • 


(42) 
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wherc  from  Eq.  (34)  and  the  section  "Calculation  of  the  Drag  from  Induced 
Velocities", 


My)  = 


r 


JO 


(^)d 


J-. 


(1  - e'2  hK)  cos  jp(y-n)|  pdp 


r * O.) 


rTT 


n 


r (nMn 


sec  W 


-2  h 


S e C 2 0 


cos  ! /.  (y  -f) ) sec  0 sin  0 d0 


J- 


oo 


(44) 


One  way  to  solve  this  integral  equation  is  by  using  a method  similar  to  that  of 
I.  Lotz2^  in  aerodynamic  wing  theory.  This  method  consists  in  assuming 
for  the  circulation  ^ ( y ) at  distance  y from  the  plane  of  symmetry  the  formula 

no 

“ (y)  = 4 b U , A2n  + j sin  (2n  + 1)  fl  , y = b cos  0 ; (45) 

n = o 

= 0 , y|  £ b 


so  that  (y)  = P(-y)  and 
oo 

r(0)  = 4bU  > (-)nA2n4j.  (45a) 

ri=o 


The  coefficients  + ^ are  to  be  determined  by  using  the  condition  (42). 

Th.  first  term  of  the  Fourier  series  in  Eq.  (45),  namely, 

I T 


r J ( y ) = 1 b U Aj  sin  t 


4 bid  A j 


(45b) 


represents  an  elliptical  distribution  of  the  circulation.  The  rest  of  the  terms 
may  be  regarded  as  a measure  of  the  deviation  between  the  actual  distribution 
and  the  elliptical. 

.substituting  Eq.  (45)  in  (43)  and  integrating  by  parts  with  respect  to  r)  , 
we  obtain  (cf.  Appendix  IV,  A) 

oo 

“(V)  = -U  V (-)n(2„*»  Aln1  , U2ntl  («)  («>) 

n=o 
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where 


rir/ 2 f'Qo 

dO 

sin  0 


-2k  t t + b/\sec  0 . / , #\\  * , w « \ , 

« "■■■  ■ 2 1'  Sin  P)  cos(t  cos0  sinO)dt, 

t - b,\scc  0 


x h J , 1 

\ = £ = depth  - -j  span  ratio  . 


(46a) 
(46b) 

In  a similar  manner,  substitution  of  (45)  in  (44)  gives  (cf.  Appendix  IV,  B) 

(JO 

'M*  •^bZ(2n*  '>  . <4’! 


Wt 


who  re 


n = o 


R>n+1(0)  sin(2n+l)0M-)n+l(3in0)R2  111*  Igi -l*2s  (2X.ljc 

\Jl  + (ZX  - icos0)2 

rv/ 2 


+ ( -)n  4b.-<  sin  0 


[• 


-2hr\scc  0 


■*2n+l^''sec*  5in  °)C3S(br^cos  dsecOsinO) 

3. 


s e 

s Ln 


c di 

;nrj dl 


(47a) 


For  symmetrical  wing  we  require  that  both  the  chord  c(y)  and  the  incidence 
Qa(y)  are  even  in  y,  hence  we  may  expand  the  following  quantities  into 
Fourier  series  with  known  coefficients 


c sin  0 

o 


GO 


TT7T 


/ £n 
n=o 


C , cos  2n  0 , 


GO 


i(y)sin0  ■ ^>^B2ncos2n0  . 


n-o 


(48) 

(49) 


Substituting  Eqs.  (45)  - (49)  into  Eq,  (42),  we  get 
- oo  - oo  .. 

rf  I XC2nC0Sin,!  zL*2ntlsin(in*1)0 

J e Ln=o  -*  41=0  J 

- l-5(-)nU"O)O2n+ll/ta+,(«'|[XBtoco.2„0] -§(^»A2ntlW2ntl(0) 

n=o  n=o  n=o  (50) 


where  ^^n+l  anc*  ^2n+l  are  8'ven  ^y  Eqs.  (46a)  and  (17a),  To  determine 
the  coefficients  ^^n+l  Eq.  (50)  would  necessitate  the  expansion  of  the 
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left  side,  ar.d  both  and  'v^  ^(0),  *ntu  a Fourier  series  in  $,  thus 

lending  to  infinitely  many  linear  equations,  in  infinitely  many  unknowns.  The 

j.lution  to  this  problem  is  obviously  n difficult  one.  However,  n very  good 

approximation  can  be  obtained  by  using  a method  for  practical  calculation  on 

28 

wings  of  finite  span,  due  to  H.  Glaue.t  . Replace  the  infinite  series  repre- 
sentation of  r (y)  given  by  Eq.  (45)  by  an  approximate  value  in  terms  of  a 
finite  series  of,  say,  m+  1 terms, 

m 

f (y)  = 4 b U 2^  A^n+j  sin(2n+l)0,  y = b cos  0 . (51) 

n = o 

This  expression  of  r'(y)  then  reduces  the  condition  of  Eq.  (42)  to 


8b 

a 


#n  r ' ■>  f rrr> 

7(77  lzL,A2n+lSin(2n+1)rf  = 1 ’E_(‘)n(2n+1)A2n+in2n+l(0) 

• Sn=o  _ - - n=o 


(aa(y)  sinff) 


m 

V 


Z.(2n+1)  ^n+l  >2n+l<*> 


(52) 


n=o 


This  equation  cannot  be  satisfied  identically  for  all  values  of  0 . However, 
if  (m  + 1)  particular  values  are  suitably  chosen  fer  0,  we  get  m + 1 linear 
equations  from  which  the  (m+  1)  coefficients,  , A^,  - - A2m+i  * can  be 
determined.  The  values  of  these  so  obtained  will  satisfy  Eq.  (52), 

not  identically,  but  only  at  the  selected  points.  In  general,  the  first  four  co- 
efficients usually  give  a sufficiently  accurate  result.  The  detail  of  such  cal- 
culations will  not  be  given  here;  however,  the  calculation  of  the  first  order 
term  is  quite  similar  to  that  of  a specific  example  of  the  direct  problem  dis- 
cussed in  the  section  "Example  - Elliptical  Distribution  of  the  Circulation 
Strength". 


Having  obtained  the  value  of  r(y)  for  this  indirect  problem,  the  cal- 
culation of  the  lift  and  the  drag  is  then  the  same  as  that  of  the  direct  problem 
discussed  in  section  "Calculation  of  the  Drag  from  Induced  Velocities"  and 
section  "Calculation  of  the  Total  Lift".  The  results  can  be  directly  written 
down  by  ascribing  to  f(^)  and  g(j.i)  the  following  value; 


P b 


- ~ 


r (q)  co s p.r)  dr) 


'-b 


sin  p r) 


d 


oo 


i^yun+l)  A2n+  I 

k?b  Jo 


n=o 

oo 


sin  (p  bcosy  )cos(  2n+  1 ) y d y 


J2n*l 


n = o 


(53a) 
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and 


g(n)  r 0 


(53b) 


MINIMUM  DRAG  FOR  GIVEN  LIFT 


>‘/e  shall  consider  a hydrofoil  of  span  2 b so  that  P'(y)  = 0 for  |yj  b. 
The  problem  is  to  find  a distribution  P(y)  such  that  the  drag  D is  minimum 
under  the  condition  of  constant  lift  L.  By  the  method  of  undetermined  multi- 
pliers the  requirement  is  ' 


• D - \ >5  L = 0 (54) 

where  the  variation  is  carried  out  such  that  ^(y)  = 0 at  y = + b.  Applying 
the  above  variation  to  the  general  expressions  for  D and  L given  by  Eqs. 
(?6a)-(29b),  (31),  (35a),  and  (35b),  using  the  relation  that  * [ f q)T (q1 ) | 

- r(q)6F(r>')  + r'(q')5P(q)  and  noting  that  both  F(q-q')  and  G(q-q')  defined 
respectively  by  Eqs.  (35b)  and  (29b)  arc  symmetric  with  respect  to  (q  -q'), 
we  obtain 


<MC,  y,  -h) 
3 z 


vo  (u  + z -hJ  )}e  r (y)  d y = 0 . 


Since  <5>  (y)  is  otherwise  arbitrary,  the  quantity  inside  the  bracket  must 
vanish  identically.  Or,  using  the  present  notation,  we  have 


2 w(y)  - \Qfu  + 2u(y)j-  0 . 


(55) 


The  constant  multiplier  \ determined  by  Eq.  (55)  equals  approximately 
twice  the  value  of  the  downwash  angle  € for  not  too  shallow  submergences; 
more  precisely,  the  difference  ( - ^o  is  a second  order  small  quantity 
given  by 


c 


u(y).  w(y. 

u(y)ifd 


+ My)i 


Therefore,  the  condition  for  the  minimum  total  drag,  accounting  for  all 
causes  and  holding  lift  constant,  is  that  the  total  downwash  angle  must  be 
constant,  along  the  span,  up  to  the  first  order  term.  When  depth  of  submer- 
gence becomes  infinite,  u(y)  tends  to  zero,  the  above  condition  is  reduced 
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to  tli o requirement  of  constant  downwash  which  is  in  agreement  with  aero- 
dynamic wing  theory. 

Substituting  the  explicit  expressions  of  the  induced  velocities  w(y)  and 
u(y)  given  by  £qs.  (43)  and  (44)  into  Eq.  (55),  we  obtain  the  following  integral 
equation  for  P (y): 


b 

p. 

r r(q)dn 

J-b 

% 

(1  - e '2h  K)  cos  p(y  -q  )>»  dp 


, TT  / 2 


+ 4 


e **  ^s*-c  ® cosl  '(y-q)  sec2©  sin  0)  sec"*©  d 0 


Jo 


- 

o b 

' oo 

>n/2  . ' 

- 

(q)d*i 

- 2h  ^ . 

e ^ u d p 

cos(p(y -q  ) sin  0)^— — — ~~T“  d C 

TT 

_ V 

-b 

Jo  J 

u - 1 sec  0 

0 J 

(55a) 


To  solve  this  equation  a practical  method  of  approximation  similar  to  that  used 
in  solving  Eq.  (42)  of  the  indirect  problem  may  also  be  applied  here.  Assum- 
ing that  the  exact  solution  of  Eq.  (55a)  may  be  approximated  by  the  expression 
given  by  Eq.  (51),  Eq.  (55)  can  be  reduced  to  the  form 


m m 

^(2n*l)  AJn<1«v2n+i(0)  --  Vinf<  jj  • ]T( -)"  U»*  1 ) 

n=o  ’ n=o 


(0) 


(56) 


where  anc*  ',n+l  arc’  8‘Vcn  by  Eqs.  (4tia)  and  (47a)  respectively. 

From  this  equation  the  firs*  ( rh+  l)  coefficients  A may  then  be  determined 
for  (m+  1)  particularly  selected  values  of  0 . It  should  be  remarked  here 
again  that  the  solution  so  obtained  cannot  satisfy  Eq.  (5b)  identically  in  P, 
bit  nevertheless  gives  a good  approximation.  The  estimation  of  the  deviation 
of  this  solution  from  the  elliptical  distribution,  which  is  the  solution  of  the 
corresponding  problem  in  aerodynamic  wing  theory,  will  be  made  in  the  speci- 
fic example  discussed  below. 
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EXAMPLE  - ELLIPTICAL  DISTRIBUTION  OF 
THE  CIRCULATION  STRENGTH 


As  shown  in  the  indirect  problem  that  the  first  term  in  the  assumed  ex- 
pansion of  P(y)  given  by  Eq.  (45)  represents  the  so-called  elliptical  dis- 
tribution, it  would  bo  of  interest  to  consider  a direct  problem  of  a hydrofoil 
of  span  2b  immersed  at  a fixed  depth  h with  .an  elliptical  distribution  of 
P(y),  that  is, 


P(  Y)  = 


I y i b . 


= o , 


I y I • 


(57) 


Although  for  a given  hydrofoil  of  fixed  geometrical  form  this  distribution 
cannot  be  held  for  different  values  of  h,  nevertheless  it  gives  interesting 
results.  In  this  problem  we  shall  further  define  the  F roude  number  of  the 
motion  by 


<r 


1 

KT 


u 

7* 


(58a) 


which  is  one  of  important  parameters  of  this  problem.  We  shall  mostly  be 
concerned  with  large  values  of  c,  corresponding  to  shallow  submergences. 
For  instance,  for  U equal  to  80  ft/sec  and  h equal  to  4 ft,  er  is  approxi- 
mately 50.  We  shall  further  denote 


P 


P is  in  general  also  very  large  in  our  velocity  range  of  interest. 
For  this  distribution  of  P(y)  we  find,  from  Eq.  (27), 


(58b) 


zr 

- -v1! 


r b 


JO 


2 JL(nb) 

\ i - cos(^y>dy = 1 0“7“ 


g(n)  = 0 


(59) 


a . Calculation  of  the  Drag 

Substituting  Eq.  (59)  into  Eq.  (26),  we  have  (cf.  Ref.  24,  p.  405): 


D, 


"pp„2r°°  Ji2(pb) 


IT  P' 


(60) 
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md  from  Eq.  (28)  we  obtain  (cf.  App.  IV,  C) 


"Pr„ 


D2  + l)3 


r 


dp 


it  p r 


-M 


I ♦ k 


Kl-L-I-  £[ 

\ 


1 + k 


+ k 


2 J 


(61) 


wlu  rc  K(k)  and  E(k)  denote  the  complete  elliptic  integral  of  the  first  and 
second  kind  respectively  end  k is  the  depth  - -j-  span  ratio  defined  by  Eq.  (58b). 
Using  the  series  expansions  end  the  asymptotic  expansions  of  K and  E 
(cf.  Ref.  26,  p.  73),  we  find  the  behavior  of  (D^  + D^)  for  small  and  large 
values  of  \ respectively  as  follows: 


D •>+  D , = 


D2+D3 
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i -i  _L_ 

tt  : g 

\i  1 + 
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as  h -*  0 ; 

as  h -» ao  , 


, 4\il+k2x 
log— — -) 


(62a) 

(62b) 


As  h-*oo,  the  drag  due  to  the  surface  effect  diminishes  like  ^ . '.hen  th. 
hydrofoil  becomes  a pi  inning  surface,  i.c.  h— *0,  we  find  that  (D,  + D,) 

l u ^ 

tends  to  a finite  value  -np  'o"/8  which  cancels  D^,  and  consequently,  the 
drag  is  the  a solely  due  to  the  wave  effect. 

To  find  the  wave  drag,  we  substitute  Eq.  (59)  into  Eq.  (29)  and  obtain 
the  following  integral*: 


°4  = ffPj“o2 


* — scc'3  i J.(i  sechstnQ) 

, »•  i . It-.  . — : 


sec  0 d 0 


sec  0 s l n W 


(63) 


From  the  integral  representation  of  the  wave  drag  given  by  Eq.  (63),  we 
note  that  ^ 

D j = p Cj  f(<r,  k),  or,  = f(<r,  k) 

p U c 

o 

which  can  be  directly  derived  from  dimensional  analysis  (cf.  Ref.  5,  p.  4 38). 
This  is  a useful  non-dimensional  form  for  experimental  purpose. 
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For  general  values  of  the  F roude  number  <r,  the  above  integral  may  be  con- 
verted into  an  infinite  series,  each  term  of  which  contains  modified  Bessel 
functions  of  the  second  kind,  (cf.  App.  IV,  D),  as  follows: 


‘V 


n-o 


For  very  small  values  of  <r  (•*  large  h,  or  X)  , we  can  use  the  asymptotic 
xpansion  of  Kn(i)  $ 


1 + 


4r. 


TT 


41  4 

o 11  /n\b 
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,2  + 


2'  (8) 


(65) 


So  we  see  that  in  this  case  the  series  in  Eq.  (64)  converges  very  fast  and  the 
final  result  may  be  expressed  asymptotically  by 


1 

8 6p2 


a + 0 (a2)  , 


as  (t-*0,  (66) 


which  diminishes  exponentially  with  increasing  h.  As  a remark,  the  above 
asymptotic  expression  of  for  a small  can  also  be  obtained  by  evaluating 

the  integral  of  Eq.  (63)  in  a complex  0 plane  and  applying  the  method  of 
steepest  descent  (e.  g.  cf.  Ref,  23,  p.  504). 


For  large  values  ol  j ( \ small  h),  Kn(— ) 6as  the  following  expan- 


sion: 
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ko(1)  = i0(i)  */  ^77T  ' <m“>  • 
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n a 2 / 


m = o 
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(67) 
(n  ^ 1) 


m = o 


It  can  be  seen  that  the  most  important  contribution  to  the  sum  of  the  series 
in  Eq,  (64)  comes  from  the  term  ^ nl  (2cr)n+*  of  the  expansion  of  ^n+j(— ). 

As  shown  in  Appendix  IV,  E,  this  term  results  in  a scries  which  converges 
only  for  Xi  1 and  represents,  in  this  region,  an  analytic  function  of  X. 
Consequently  the  value  of  the  wave  drag  for  shallow  submergences 


t 
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(~  \ ' 1 ) may  be  obtained  by  some  analytical  continuation  of  this  function  from 
the  region  \ 1 to  the  whole  region  of  physical  interest,  namely,  X.  '•  cd 

(cf,  App.  IV,  E).  The  final  result  is 


t" ~t~ I i-hJ 
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(68) 


It  is  of  interest  to  note  that  when  the  hydrofoil  is  close  to  the  surface,  the 
wave  drag  is  rather  similar  in  behavior  but  of  the  opposite  sign  and  numeri- 
cally almost  twice  of  the  drag  due  to  the  surface  effect  (cf.  Eq.  (61)).  How- 
ever, it  should  be  emphasised  here  that  even  1 * (y)  were  the  same  function 
of  y for  all  values  of  X,  l“  still  depend  on  X,  especially  for  X small, 

(cf.  Eqs.  (9)  and  (41a)).  Consequently  the  actual  dependence  of  the  induced 
and  wave  drag  on  X for  a given  hydrofoil  at  fixed  geometric  incidence  can 
only  be  determined  after  the  function  rQ(X)  is  found  later  in  this  section, 
part  (d). 


It  is  also  of  interest  to  compare,  when  a is  small,  the  first  term  in 
Eq.  (64), 

it  p P 2,2  « * , . 

et"  ! Ko<;>  * (l*f)  Ki<?> 


TT 


with  the  wave  resistance  on  a moving  sphere  of  radius  r given  approximately 
by  (cf.  Ref.  5,  p.  437,  also  Ref.  6): 
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( - U \ 
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Except  for  a proportionality  constant  which  depends  on  the  size  of  the  obstacle, 
they  have  the  same  dependence  on  <r.  Evidently  an  observer  on  water  sur- 
face cannot  distinguish  bodies  of  different  forms  if  they  move  at  sufficiently 
great  depth. 
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b.  Calculation  of  the  Lift 

Fho  tot.il  lift  experienced  by  the  hydrofoil  can  be  obtained,  by  substitu- 
ting Eq.  ('9)  into  Eqs.  (M)-{35),  as  follows: 


L = L + AL.  + Al., 
o 1 2 


who  re 


L r p U r 
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(72) 


It  c~n  be  seen  that  L(  is  the  aerodynamic  value  of  the  lift,  AL  j is  inde- 
pendent of  the  gravitational  effect  and  thus  represents  the  correction  due  to 
the  mean  free  surf'ce.  All  the  lift  caused  by  the  surface  wave  effects  is 
represented  by  aL^.  For  small  incidence  ingles,  is  approximately 

proportional  to  the  incidence  angle  (cf,  Eq.  (4  1a)),  hence  relative  to  Lq  , 

Al.  and  D are  second  order  small  quantities. 

The  integral  representing  AL  ^ contains  the  same  integral  given  in 
Eq.  (hi)  (cf.  also  App.  IV,  C),  the  result  of  which  may  be  then  applied  here. 
AfL  r this  substitution,  integrating  by  parts  with  respect  to  0,  we  obtain 
: simpler  representation  for  AL  j (cf.  App.  IV,  F)  : 
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wh.rc  w(k)  is  - derived  complete  elliptic  integral  (cf.  Ref.  2S,  p.  73) 
de  fine  d by 
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and 


-X 


(73b) 


M 1 + 


For  0 < h v oo  we  have  the  corres  ponding  value  s of  k and  > as  0 < k < oo 
and  1 >>>  0.  When  k— * 0 (X- *1),  integral  in  Eq.  (73)  tends  to  a definite 
value,  namely 
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where  B(k)  is  anothe  r de  rived  complete  elliptic  integral  (cf.  Ref.  25,  p.  78) 
defined  by 


B(k)  = [ E(k)  - (1-k2)  K(k) 


(74a) 


Hence  Eq.  (73)  indicates  that  aL.  tends  to  CD  like  i as  \— *0,  noting  that 
Z l \ 

i*'  docs  not  vanish,  as  will  be  shown  later.  This  fact  is  not  surprising 

because  it  is  known  that  as  the  hydrofoil  approaches  the  planning  condition, 

the  total  lift  drops  to  almost  half  of  its  aerodynamic  value  Lq  (e.g. 

cf.  Ref.  14),  The  divergence  of  the  integral  in  £q.  (71)  simply  implies 

that  this  second  order  small  quantity  will  grow  so  large  that  it  will  modify 

the  value  of  the  first  order  quantity  Lq. 

For  X.  not  too  small  (~x  not  lo°  close  to  1),  the  integral  in  Zq.  (73) 
can  be  calculated  by  using  the  known  expansion  of  *(k),  (cf.  Ref.  25,  p.  73), 
given  as  follows: 
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Substituting  this  expansion  into  Zq.  (73),  then  applying  the  transformation 

2 2 

k =>  t and  integrating  termwise,  we  obtain  (cf.  App.  IV,  G)  : 
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This  formula  should  be  good  for  X,  > 1 (or  h>b).  For  f}  close  to  1 (h£;0), 
the  above  series  converges  slowly;  hence  we  have  to  resort  to  some  other 
approximation  in  order  to  facilitate  calculation.  For  this  case,  we  define 
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Now  the  first  integral  can  be  integiatcd, 

KX 

c$(k)dk  = j[>.B(x)  -XX'2  CSJX)]  . 
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and  the  second  integral  may  be  appioximated  by 
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Finally  we  have 
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and  using  the  known  expansion  of  BpO  and  (S(X)  as  we  obtain 
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+ 0(X  log  X)  as  X-* 0 (78a) 


The  integral  representing  in  Eq*  (72)  is  indeterminate,  but  its 

principle  value  exists.  The  method  of  contour  integration  described  pre- 
viously >n  the  section  "General  Formulation  of  the  Problem"  gives  a complicated 
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expression.  In  this  case  we  shall  use  the  following  approximation.  By 
p = (<sec  9)u  we  can  transform  Eq.  (72)  into  the  following  form: 
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The  inner  integral  is  very  similar  to  that  representing  in  Eq.  (63)  which 

is  discussed  in  detail  in  Appendix  IV,  D.  This  integral  may  be  treated  in  a 
similar  way  (cf.  App.  IV,  il).  Because  the  interesting  feature  of  aL^  also 
comes  from  the  integration  with  respect  to  u , wo  shall  take  only  the  first 
order  term  of  the  .nner  integral.  The  final  result  is  as  follows(cf.  App.  IV,  H) 

for  c very  large  (shallow  submergence), 
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and  for  o very  small, 
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From  these  relations  it  is  of  interest  to  see  that  AL-,  tends  to  a constant 

2 j “■ 

value  p'*’  P (— )/-i\Z  as  h tends  to  iero.  This  result  is  quite  a contrast  to 

the  feature  of  ALp  Compared  with  AL  ^ , the  value  of  AL.^  is  smaller  for 
o large,  and  greater  fur  o small. 


c . The  Induced  Velocity  at  the  Hydrofoil 

In  order  to  investigate  further  the  geometric  and  hydrodynamic  proper- 
ties of  the  hydrofoil,  we  need  to  know  the  value  of  Pq  ; and  if  we  want  to  cal- 
culate in  terms  of  given  quantities,  we  have  first  to  obtain  the  values  of 

the  u-velocity  .and  the  downwash  at  the  hydrofoil  (cf.  Eq.  11). 

Substituting  the  value  of  P(y)  into  Eq.  (4  3),  we  may  express  u(y)  in 
terms  of  two  parts  : 

n(y)  = u^y)  + “2(y) 


(82) 
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The  integral  in  Eq.  (82:.),  which  represents  the  effect  due  to  the  mean 
free  surf-tee,  cm  be  evaluated  exactly.  The  result  is  a s follows, 

(cf.  App.  IV,  1)  ; 

00  n 

r o 


^(y)  = - 


9 jT  b 


1 \^(-)’(n+2)  2n  11  1 

__  > F(n+  n;2; ^ ) 

. '•  (1+1  \ ) 1 + ix 


n b X ; \‘  n = o 


(83) 


The  v-’lue  of  Uj(y)  is  •>.lwnys  negative  and  becomes  infinite  ns  h-»0.  It  can 
be  ihjwr,  that  the  above  series  converges  uniformly  with  respect  to  q when 
OS  ri  c.  1 V/e  inav  also  find  that  j ( V ) is  a slowly  varying  function  of  y ; 

m particular,  the  value  of  Uj(0)  is 
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where  B(k)  is  a derived  complete  elliptic  integral  defined  by  Eq.  (74a). 

11  can  be  ->howr.  tu  t u,(y)  is  also  a slowly  varying  function  of  y,  i.e. 
the  value.'  of  u (y)  at  points  c lo  e to  the  hydrofoil  differ  very  slightly  from 
its  va.lue  at  plane  o;  symmetry  t>iVen  by  (cf.  App.  IV,  J)  : 
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and 
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It  is  of  interest  to  see  that  Uj(0)  and  uz(0)  are  of  the  same  order  of 
magnitude.  Therefore  the  totai  u-component  velocity  at  y = 0 is 
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The  downwash  at  the  hydrofoil  can  be  calculated  in  a similar  way: 
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and 
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Now  the  integral  in  Eq.  (87c)  can  be  evaluated  while  that  in  Eq.  ( 8 7d ) again 
has  to  be  approximated. 
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It  can  be  seen  that  Fj(X,rj)  is  a slowly  varying  function  of  q for  fixed  \, 
In  particular,  we  have 


Fj(X,  0)  - 


i|  2 X 


'll  1 + 4X2 


(89a) 


F,(X.  + 1)  = 


X \|  X 'li  + x2  + x2  + 1 x Lv2  - x2 


(89b) 


X \ 1 + X‘ 


For  small  values  of  a,  application  of  the  method  cf  steepest  descent  to  the 
integral  in  Eq.  ( 8 7d ) yields 
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which  diminishes  exponentially  with  increasing  h and  is  independent  of  y for 
points  close  to  the  hydrofoil.  For  very  large  <r,  the  integral  in  Eq.  ( 8 7d ) also 
represents  a slowly  varying  function  of  y and  consequently  only  the  value  of 
<p„  (0,0,  -h)  is  of  interest.  It  can  be  shown  that  (cf.  App.  IV,  K): 

tz 


(0,  0,  -h)*- 
•z 


rQ 

71 


l* 


“V  - 


2 4 ( *2  * 


(0,0,  -h) 


(91) 


Finally  we  have  the  value  of  w(0),  by  substituting  Eqs.  (87a),  (88)-(91)  into 
Eq.  (87a),  as  follows 
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which  tends  to  aerodynamic  value 
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as  h — * oo  . 

It  may  be  remarked  here  that  both  u(y)  and  w(y)  are  almost  constant 
spa.nwise  at  points  close  to  the  hydrofoil.  This  fact  justifies  some  of  the 
basic  assumptions  introduced  previously  in  the  sections,  "General  Formulation 
of  the  Problem"  and  "Geometry  of  the  Hydrofoil;  Effective  Angle  of  Attack". 
This  result  also  supports  the  .assumption  that  the  same  circulation  distribution 
will  hold  for  a wide  range  of  depths. 


d.  Geometry  of  the  Hydrofoil 

We  are  interested  primarily  in  the  case  of  shallow  immersion,  because 
the  geometry  of  the  hydrofoil  at  deep  immersion  will  be  the  same  as  the 
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cor  re  sponding  aerodynamic  problem.  Substituting  Eqs.  (86)  and  (92)  into 
(4  1a)  and  neglecting  aecond  order  small  quantities,  we  obtain  the  value  of 
Pq  (:P(0)  ) a a a function  of  X.  and  c as  follows: 


P 

o 


1 4- 


1 


c 

o 


a 

c 


f2(\, a ) ♦ 


2 

TT 


a f 

a 


(93a) 


w here 


f j(^.  <r)  - 


and 


2Bf 


b (-drrr)  ♦ 


V tt 


4\‘{7+4?  \ I + -! \2  ^2K(U4\Z)1/4  {lVj\2+2v] 


+ °v 

(93b) 


f2<^)=  [l  - j=]  l,+0<7* 

1 1 4 \ 


(93c) 


/s  h tends  to  oo,  f^  tends  to  zero  and  f^  tends  to  ; hence  the  value  of 
at  infinite  depth  becomes 
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which  is  the  well-known  aerodynamic  value  of  I-1  in  terms  of  given  quantities. 
It  is  of  interest  to  note  that  changes  with  respect  to  P (0),  as  X.  and 

a vary,  according  to  the  following  relation: 
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which  is  a function  of  \,  <r  and  a for  a given  hydrofoil.  This  function  and 

2 3 

its  squared  value,  y-  , are  plotted  against  \ in  Fig.  3 for  two  incidences 
Qa  = 3°  and  6°  under  the  following  operating  condition- 
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The  values  of  y and  y*  decrease  with  decreasing  and  tend  to  unity 
asymptotically  as  \-*oo  . For  small  values  of  X.  (\  < 1),  both  and  Pq2 

deviate  appreciably  from  their  respective  values  at  infinite  depth. 


It  may  be  pointed  out  here  that  because  both  u(y)  and  w(y)  are  almost 
constant  along  the  span  of  the  hydrofoil,  it  follows  from  Eq.  (4  1)  that  for  an 
elliptical  distribution  of  P(y),  the  plan  form  is  also  nearly  elliptical, 

i*4-  <97) 

b 

together  with  a negligible  geometric  twist. 

Another  consequence  of  the  above  results  is  that  the  condition  of  Eq.  (55) 
for  minimum  drag  with  prescribed  lift  is  also  approximately  satisfied  by  almost 
constant  (spanwise)  u(y)  and  w(y).  Hence  the  hydrofoil  of  elliptical  plan  form 
still  gives  approximately  the  minimum  drag  for  a given  lift. 


c(y)  s Co  \ 


e . Over-all  Hydrodynamic  Properties;  Lift  and  Drag  Coefficient 

In  this  section  we  shall  again  only  consider  the  case  of  shallow  submer- 
gence. Combining  our  previous  results  we  obtain  the  total  lift  and  drag  of 
the  hydrofoil  at  small  depths  as  follows: 
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and  ( 


and  ^”(\,j,a  ) are  given  by  Eqs.  (94)  and  (95)  respectively. 
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(99a) 


(99b) 


At  this  stage  we  should  be  able  to  see  the  properties  of  the  total  drag 

exerted  on  the  hydrofoil.  As  mentioned  before,  when  X — ♦ oo  D tends  to 

the  aerodynamic  value  D ■ tt  p P /8  , hence  we  write 
7 CO  y o' 
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oo 

This  ratio  is  plotted  in  Fig.  4 for  o =6°  and  under  the  operating  condition 
given  by  Eq.  (96a).  The  first  term  on  the  right  hand  side  represents  the 
contribution  from  the  total  induced  drag;  this  quantity  tends  to  asymp- 

totically for  large  h and  decreases  rapidly  as  the  hydrofoil  approaches  the 
surface  (h  < ?.b).  The  second  term  represents  approximately  the  wave  drag 
which  is  negligibly  small  for  h > 2b  and  becomes  important  only  at  small 

depths.  As  a whole,  the  total  drag  D deviates  only  slightly  from  D for 

b 1 ® 

h y ^ , and  decreases  mode  r ately  at  small  depths  (h  <.  ^ b)  which  is  ac  tually 

the  range  of  interest.  .Vhcn  the  hydrofoil  approaches  the  planning  condition, 

the  total  dr^g  is  solely  due  to  the  wave  effect. 

In  order  to  meas  ire  the  total  drag,  including  the  wave  drag,  exerted  on 
a hydrofoil,  the  experiment  must  be  carried  out  in  a towing  tank.  Then  it 
would  be  convenient  to  define  the  over-all  lift  and  drag  coefficient  with  re- 
spect to  the  towing  speed,  such  as, 


C 


L 


T 

7 


L. 


r 

WD 


D 

|7urI 


(ICO) 


where  p is  the  density  of  the  liquid,  U the  towing  speed,  and  5 the  area 
of  the  plan  form.  The  above  definition  for  and  were  used  in  pre- 

vious experiments  (known  to  the  author).  In  order  to  compare  with  these 
available  experimental  results,  we  shall  also  adopt  the  above  definition  of 
Cj^  and  Cp.  Substituting  Eqs.  (94),  (98)  and  (99)  into  (100),  we  obtain 
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cL  = J ^ (-£) 


a a 

cap 

(l  + IT  ae> 


ST  * «V 


a a 
e a f 

, . i Co 

1 + IT  V 


{fj<M-f4(*.<r)J  , 


(101) 


:d  = £ * Hr)' 


a q y- 
e a } 

1 + IT  ac 


|f5(X)  ♦ 2 [l  - f5(X)]|  . 


(102) 


where  /R  = •— t—  = aspect  ratio  and  t^,  f^,  are  given  in  Eqs.  (98)  and  (99). 
The  lift  drag  ratio  is 
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, c 
1 t °\ 

t (nr) 


a q V 

£_jlL_ 

. i c° 

1 * 5T  \ 


|£3(k)  - f4(x*  <0  j 
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. c 

T(T} 


a a v^- 
e a ; 


u + 


1 Co 

IT 
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(103) 


f5(X)  + £‘l-f5(X)_' 


U 

These  results  should  be  good  for  K (2  (h  < 2b)  provided  that  a - r-  is  large, 

a n 

say,  greater  than  20.  From  the  expressions  of  L and  D for  h large,  we 
note  that  as  h-*oo,  both  and  tend  asymptotically  to  their  aerodynamic 

values,  namely, 


;l*j  ^Hr>  ".V 

. 1 o 

1 + I T ae 


7-T  * (^)2 


'D-  FT 


a a 


• i C° 

UITa, 


so  that 


cn~—  C.  2 

D tt  *R  L 


which  agrees  with  aerodynamic  wing  theory.  However,  it  can  also  be  seen 
from  Eqs.  (101)  and  (102)  that  even  for  \ small,  C^/C^  *s  stHl  propor- 
tional to  1/  /‘R  if  the  small  quantity  of  the  second  order  is  neglected.  From 
this  result  it  follows  that  in  order  !o  improve  the  hydrodynamic  properties  the 
preference  should  be  for  a high  aspect  ratio  hydrofoil,  even  for  operations 
near  the  water  surface. 

To  illustrate  the  details  of  the  behavior,  and  to  compare  our  results 
with  experiments,  we  shall  take  specific  values  of  cq,  b,  >R , and  ag  to 
determine  all  the  coefficients.  This  will  be  shown  as  follows. 


«>  . 
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DISCU5SION  OF  RESULTS 


Consider  a hydrofoil  of  elliptical  plan  form  with  other  specifications 
and  an  operating  condition  given  by  Eq.  (96a)  so  that  - 6.  3.  These  data 

serve  to  determine  the  coefficients  in  Eqs.  (101) -(103)  which  then  become 


CL  =4.74c 


af  °*  306  aaV2  {f3(V)  ' f4(K'a)}  * 


CD  = 1.  14  aaZr2  {f5(K)  4 2[l  - f5(X)]} 


(104) 


(105) 


1 - 0.065  aar{f3<M  - UV.cr)  } 

°.^Qar  ff5(M^  [ l-f5(\)j] 


(106) 


These  equations  '.re  plotted  against  X.  for  = b in  Fig.  5 and  also 
plotted  against  for  several  small  values  of  \ in  Fig.  6.  From  these 

curves  several  interesting  conclusions  may  be  drawn. 

/ t depths  greater  than  2 chords,  the  influence  of  the  surface  of  the 
water  is  negligib'.y  small  and  the  hydrofoil  will  have  characteristics  similar 
to  those  of  an  airfoil  of  the  same  section.  In  the  range  of  depths  less  than 
2 chords,  Cp  decreases  gradually  oecause  the  decrease  of  liquid  flow 
(with  respect  to  U)  auove  the  foil  diminishes  the  induced  dovnwash.  At  very 
small  depths  (h  < cq/2),  the  wave  drag  only  is  important,  our  results  show 
that  when  the  hydrofoil  is  ne..r  the  surface,  the  rate  of  energy  shed  to  form 
the  wave  system  .s  slightly  less  than  the  induced  drag  at  large  depths.  In 
the  same  range  of  operation  (h  < 2 c ),  C.  decreases  comparatively 
rapidly  to  almost  C.  as  the  hydrofoil  approaches  the  surface.  This  re- 

C L-oo 

duction  in  lift  results  from  the  decrease  of  mass  of  water  flowing  over  the 

upper  surface  of  the  hydrofoil,  causing  a reduction  of  the  absolute  value  of 

the  negative  pressure  on  the  suction  side.  The  corresponding  lift-drag 

ratio  decreases  very  slowly  with  decrease  in  depths  for  h > c / 2 , and  is 

almost  constant  about  the  point  h - -y  cq.  This  ratio  decreases  rapidly  with 

further  decrease  in  cLpth  for  h < -L  c . This  result  also  indicates  that  if  the 

7 ° l 

water  surface  is  not  too  choppy,  then  the  depths  between  cq  and  1 cq 
would  be  a favorable  range  for  operation.  This  optimum  range  of  depths 
also  associates  with  a stabilising  effect  because  this  range  corresponds  to 
the  middle  part  of  increasing  slope  dC^  /d\  , consequently,  any  further 
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decrease  in  depths  will  result  in  a rapid  decrease  of  lift  so  the  hydrofoil  will 
sink,  and  if  it  sinks,  an  increase  in  lift  will  raise  it  up  again.  The  above  re- 
sults as  shown  in  Fig.  5 agree  well  with  some  of  previous  experiments  (cf. 

Ref.  >7,  18,  4). 

If  we  add  an  almost  constant  frictional  drag  (about  = 0.02)  to  the 

drag  coefficient,  then  the  result  shown  in  Fig.  6 is  m good  agreement  with  the 
observations  made  in  Ref.  21  and  18.  The  reduction  in  and  Cp  with  de - 
c rease  in  depths  becomes  more  appreciable  at  larger  angles  of  attack. 

It  is  also  of  interest  to  consider  the  situation  of  practical  operation. 

When  a hydrofoil  moves  through  water  of  infinite  surface  extent  at  shallow  sub- 
mergence, an  observer  who  moves  with  the  hydrofoil  can  measure  the  value  of 
the  approaching  flow  velocity  U ♦ u{0)  more  easily  (for  instance,  with  a pitot 
tub^.)  than  the  value  of  U at  upstream  infinity.  Then  if  one  defines  the  over- 
all lift  and  drag  coefficient  based  on  the  value  of  U + u(0),  that  is 

c ’ = 2 (107) 

u 1 .,2.,  A u(0).  2 . 

p U ( 1 + — 'jj')  o 


|rv*.<r)]  • {108) 

one  will  find  that  both  CjJ  and  C^'  are  almost  constant  for  X > 0.  1 (cf. 
Fig.  7)  although  Cjj/C^'  is  identical  to  C^/C^.  However,  the  total  lift 
L and  the  total  drag  D calculated  from  Eq.  (107)  would  still  have  the  same 
dependence  on  X as  and  C^,  only  with  a different  proportionality 

constant. 

The  above  results  are  derived  under  the  assumption  of  large  Froude 
numbers,  the  effect  of  speed  U on  and  is  dropped  out  from  our 

final  formulas.  This  effect  becomes  significant  for  motions  with  small 
Froude  number  (say,  less  than  1),  For  instance,  let  us  consider  the  same 
hydrofoil  as  that  in  the  previous  case,  but  moving  at  a low  speed, 


C ' = 


7 pU2(l  + ^)2S 


whe  re 


1 

77, 


(-£-)  a o 

p e a 

l 1 C° 

1 + 5 "F  3e 


c = 8ft, 
o 


a = 2 tr 
e 


U = 10  ft/sec, 


b = 20  ft, 


(109a) 
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then 


0.  156  , 


0.  156 

«r  = — y— 


which  may  be  considered  as  small  for  X > . Using  the  values  of  u(0)  and 

w(0)  for  <r  small  (cf.  Eqs.  86,  92  and  41a),  we  obtain  the  value  of 
f{\,  <r,  a ) in  this  range  of  operation  to  be 


f(X.  <r,  oa) 


where 


1 + 


1 Co 

¥ T 


a 

e 
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I 

j 


(o'  <\) 


(110a) 


(110b) 


g2(X)  = 2 


X 

U+4\2 


(110c) 


The  value  of  y given  by  Eq.  (110)  is  plotted  in  Fig.  8 for  = 6°,  using 
the  specification  given  by  Eq.  (109).  The  resultshows  that  f'  also  decreases 
with  decreasing  X for  «r  small,  but  the  rate  of  decrease  is  slower  than 
that  in  large  <r  case. 


Since  the  F roude  number  effect  only  enters  in  the  expression  of 

and  AL^.  it  suffices  to  discuss  only  the  behavior  of  these  two  quantities. 

With  respect  to  a convenient  reference,  say,  D , we  have  the  values  of 

D.  and  aL,  for  a small  as  follows: 
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(112) 


Noting  that  p is  proportional  to  U , it  can  be  seen  from  above  relations 
that  for  small  Froude  numbers  (-v  low  speeds)  an  increase  in  speed  will 
cause  both  and  to  decrease,  with  the  drag  decreasing  more 
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rapidly.  This  effect  has  also  been  found  experimentally  at  small  values  of  U 
(cf.  Refs.  17,  21).  £qs.  (Ill)  and  (112)  are  plotted  in  Fig.  9 with  the  spccifi 
cation  given  by  (109).  In  the  present  range  of  operation  (it  = — , K >*j) 

the  wave  drag  is  negligibly  small  compared  with  , but  aL ^ increases 

very  rapidly,  relative  to  D^,  with  decrease  in  depth.  The  total  change  in 
in  this  range  is,  however,  still  quite  small  with  respect  to  Lq. 
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Fig.  7 - Effect  of  depth  on  C^'.  C^'  , 

(based  on  approaching  flow  velocity). 


Fig.  8 - Effect  of  depth  at  small  values  of  Froude  number. 
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appendi:;  i 

Some  Integral  Representations  of  a Lifting  Line 

It  is  given  that  a lifting  line  of  span  2b  is  located  along  y-axis  from 
-b  to  b with  a known  distribution  of  circulation  |~‘ ( y ) , (see  Fig.  9),  the  free 
stream  being  uniform,  of  velocity  U in  x-direction.  The  problem  is  to  find 
the  induced  velocity  potential  <p(x,y,  z)  due  to  this  lifting  line. 


Fig.  9 - A lifting  line  with  its  trailing  vortex  sheet. 

First  this  potential  y should  satisfy  the  Laplace  equation 

VZ*  = 0 . (I.  1) 

Next  wc  shall  impose  boundary  conditions  for  this  p.oblem.  According  to 
the  lifting  line  theory,  the  approximation  may  be  made  that  the  trailing 
vortex  sheet  is  parallel  to  the  free  stream;  that  is,  on  the  surface  z = 0, 
with  x > 0,  jyj  ,<  b.  It  is  easy  to  see  that  boundary  conditions  for  <p  should 
be  as  follows  (see,  for  example,  Ref.  22)  : 
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(i)  <P  = 0 on  z = 0,  outside  the  vortex  sheet; 
(ii)  Across  the  vortex  sheet,  <p  has  a jump  P(y) 
(x,  y,  0 + ) = + Ofl  on  the  vortex  sheet; 


(1.2) 


(iii)  arc  continuous  and  discontinuous  across  the  vortex  sheet, 

' 3z  dx  gy 

but 


|d£ 

|dy 

+ 

o 

II 

3y 

z - 0 • 


, so  that  the  pressure  is  continuous  there. 


(iv)  <p  and  grad  ip  -*  U as  '^y^+  z^  — * + oo  and/or  x — * - oo  . 


Conditions  (i)  - (iv)  also  imply  that 
('>•)  <?  (x,  y,  z)  --  -0  (x,  y.  -z)  ; 

(v»)  iz  <p (x-  y* z)  = * (*.  y.  *z)  ; 

pi  V \ 

(vii)  (p(x,  y,  + 0)  = + (1  + sign  x)  . 


One  way  to  solve  (I.  1)  together  with  the  boundary  conditions  (I.  2)  is  by 
a Fourier  transformation.  De.'inc  the  double  Fourier  transformation  of 
<p  (x,  y,  z)  with  respect  to  n and  y by 


If  (K,  n,  z)  = 


r»oo 


-oo 


- i \ x 


dX 


P CD 


e * 1 ^ Y V (x,  y , z)  d y 


-oo 


so  that  the  inversion  formula  is  given  by 


<p(x,  y,  z)  = y- 


■oo 


-oo 


r oo 


i Kx  . 
e d X 


'V 

el  K y <f  ii,  z)  dn 


•oo 


(I.  3) 


(1.4) 


Applying  (I.  3)  to  (I.  1),  we  have 


,£  ,2  2.  ^ n 
( — ^ ^ n - o • 

d z 


The  solution  which  satisfies  (iv)  and  (v)  has  the  following  form 


•6  > 
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\ , A,.  v - i X2+p2  |z| 

<p  (X,  p,  z)  = (sign  z)  A(X,  p)  e r 

so  that 

<p  (X,  p,  0 + ) = A(X,  p)  . 

Application  of  the  same  transformation  (I.  })  to  (vii)  gives: 


' oo 

■>  oo 

' oo 

<P  (X,  p,  0 + ) = 

(t))*'4^  dr, 

-iXx  . 

e sign  x dx  + 

-iXx  , 
e d x 

V. 

-oo 

l-oo  w 

1-00 

Using  the  conventional  method  of  summability,  we  have 


’'CO 

-iXx  , 2 , 

e sign  x dx  = — lim 

1 i -»0  + 

- OO  sj 


■1A 


e x sin  X x d x = 


_2_ 
l X 


and 


I 


oo 


e‘Ux  dx  = 2 ttS(X)  . 


'-OO 

Hence  we  also  have 


<p  (x.p.o+)  = -L  | ^-L  + tt  <5 (x ) J 


r oo 


i (r;)  e*  1^r)  d-! 


- oo 


Comparison  of  (I.  7)  with  (I.  6)  gives 


a(x'^=T?[tx  + tt6(x)] 


1 oo 


r(n)  e'1^  dn 


-oo 


Substituting  (1.8)  into  (1.5)  and  then  applying  the  inversion  formula  (I. 
obtain  the  solution 


9 = sign  z 


r>  oo 


1 

2/ 


r oo 


f'fajd' 


-oo 


r oo 


sin  X x 
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, v - \l  xV  | z|  . 

coSji(y-Ti)  e * r 1 ' d(i 


r>00 


4 IT 
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■oo 


'OO 


cos  p (y-r) ) e ^ ,Z|  d | 


(1.5) 

(1.6) 


(1.7) 

(1.8) 

4),  we 

(1.9) 


This  form  of  integral  representation  for  is  the  one  given  by  von  Karman 
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(Ref.  22).  A translation  of  the  origin  in  the  z-direction  then  gives  Eq.  (10) 
of  the  text.  Another  convenient  form  of  (p  may  be  obtained  by  introducing 
the  following  transformation  of  variables: 

\ = k cos  9 , ^ = k sin  0,  (I. 10) 

If  one  uses  the  relations: 


r-n/Z 


cos(kx  cos  0)  cos(ky  sin©)  d©  , 


do 


and 


poo 


cos  k y = k 


Jo  (k  ^2*y2)  d K . 


yjo 


then  Eq.  (1.9)  is  converted  to: 


(I.  11) 


(1.  12) 
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or, 
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T 4 TT 
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<JO  UO 


c k ,7'1  Jo(k  ^ (x-£  )2+  (y-  rj2)  kdk 


(I.  M) 


The  choice  of  (1.9)  or  (1.  14)  is  only  a matter  of  convenience.  A close  form  of 
<p  may  be  deduced  from  (1.  14)  by  noting  that 
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(I.  15) 


(I.  16) 
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It  may  be  remarked  that  these  results  may  also  be  obtained  from  the 
differential  point  of  view  by  considering  the  velocity  field  induced  by  an 
element  of  circulation  T(y)dy  at  y and  by  integrating  it  along  the  lifting 
line. 

If  | grad  | <<  U and  the  effect  of  gravity  may  be  neglected,  then  the 
linearized  pressure  is  given  by 

p - pQ  = - p U <?x  (I.  17) 

where  from  (I.  16) 

».■  * 


p,„) 

- m it+fv-nl  + 7 I 


(1.18) 
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APPEND1X  II 

Formulation  of  the  Nonstationary  Flow  Problem  and  the 
Jtationary  Flow  as  a Limiting  Case 


Take  axes  Ox  and  Oy  in  the  undisturbed  surface  of  deep  water,  and  Oz 
vertically  upward,  the  origin  being  at  rest  with  respect  to  the  fluid  at  infinity. 
Consider  now  the  flow  caused  by  a wing  which  starts  to  move  at  t = 0 from  the 
position  (0,  0,  -h)  in  the  negative  x-direction  with  velocity  U.  The  effect  of 
an  element  of  wing  span  A r]  at  y = rj  of  circulation  f'(q)  Ary  is  to  produce  an 
initial  impulse  symmetrical  about  (0,  rj,  0)  on  the  initially  still  surface.  These 
conditions  correspond  to  initial  data  for  <p(x,  y,  z,  t)  and  4 (x  , y,  t)  as  follows: 


' p A0  (x,  y,  0,  0)  = aF  (to)  , t)2  = x“  + (>  -*1  )2 

. A 4 (x,  y,  0)  = 0 , 


(II.  1) 


where  aF  (a)  should  be  equal  to  the  initial  pressure  perturbation  on  the  free 
surface.  Using  (I.  17)  and  (I.  18)  and  considering  the  effects  due  to  both  the 
wing  and  its  image,  we  obtain 


a f (w)  = - 


= . Jl , (II. 2) 

V + (yV  + h2]}/2  ^ [to2  + h2]3/2 


The  solution  for  <p  which  satisfies  the  Laplace  equation  with  initial  condi- 
tions (II.  1)  is  (Ref.  5,  p.  431): 


r oo 
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AF  (a)  JQ(k  a)  ado  , (z  $ 0) 


and 
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\Jgk  sin(  \Jgk  t)J  (ku)kdk 


A F(a)  JQ(ko)  ado  . 
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Howc  ver , 
poo 


aF(q)  J^(ka)  ada  = 


f (r!  ) d*! 


JO 


■oo 


q J (k  a) 


d q 


£girr>',')AT1  fir  K i (kh)  (Ref-  26-  p*  30) 
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e ^ sin(  \[i"k  l)  JQ(kui)  kdk 


The  resulting  values  of  0 and  C due  to  the  entire  wing  may  then  be 
obtained  by  integrating  A0  and  AC  along  the  span: 
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(11.  3) 


It  should  be  pointed  out  here  that  these  are  only  the  results  due  to  an  impul- 
sive motion  of  a hydrofoil.  To  extend  them  to  the  case  of  a continuously 
moving  wing  in  ( -x)-direction  with  velocity  U,  we  have  to  superimpose  all 
instantaneous  disturbances  at  x = - Ut  at  time  t = t.  We  replace  t by 
t -T  and  x by  x - £ = x + Ut;  and  integrate  with  respect  to  t over  the  time 
during  which  the  system  has  been  in  motion. 


poo 


Pt 


pCO 


1 

7 ~ 


PhJdT) 


U(t) dT 


ek(z-chJsjgk(t-T)  J (k  ^(x+Ut)2*  (y-^n  kdk 


1-00 


Jo 


Jo 


(11.4) 


-56- 


4 = - 
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^jgk  sin  ^gk(t-T)JQ(k  ^ (x+U-r^+(y-r)  f)  k d k 

(II.  5) 


These  are  the  integral  representations  of  the  solution  of  nonstationary  motion 
of  a hydrofoil.  This  0 is  only  the  part  of  the  perturbation  potential  due  to 
effects  of  free  surface  and  wave  formations,  not  including  the  contribution  of 
the  wing  and  its  image. 


Now  we  shall  pass  on  to  the  limiting  case  such  that  the  wing  has  been  in 
motion  since  X = - ou  with  constant  velocity  U.  A difficulty  arises,  however, 
because  the  integral  with  respect  to  a from  -oo  to  t is  nonconve rgent.  Since 
the  experimental  evidence  ensures  that  £ should  be  finite,  we  must  seek  a 
remedy  in  order  to  agree  with  observation.  Mathematically  it  is  usual  in 
these  problems  to  introduce  a convergence  factor  e ^ and  then  let  p 
tend  to  zero  afterwards. 
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(II.  6) 
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(II.  T) 


Physically  it  is  usually  explained  that  the  indeterminateness  of  the  problem  is 
due  to  the  absence  of  viscous  effect.  To  avoid  this  difficulty,  Rayleigh  sug- 
gested the  assumption  that  the  fluid  element  is  also  subject  to  a frictional 
resisting  force  proportional  to  the  relative  velocity.  Then,  instead  of  £q.  (4) 
we  have  (Ref.  5,  p.  399) 


-E-+  -jjj*-  + gz  - P0  = 0 (II. 8) 

From  this  it  can  be  shown  that  the  circulation,  hence  also  0,  has  a factor 
c ^ l,  which  shows  the  damping  due  to  viscous  effect. 

It  is  natural  to  expect  that  expressions  (11.  6)  and  (II.  7)  will  become  in- 
iependent  of  t with  respect  to  a coordinate  system  moving  with  the  hydrofoil. 
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Hencc  we  introduce  the  following  Galilean  transformation: 


x = x + U t , 


y = y 


Z = Z , 


t = t 


(II.  9) 


and  change  the  variable  of  integration  from  to  T = t - T so  that 
x + Ut  = x - U T 

and  then  dropping  the  bar,  we  obtain: 
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J-oo  o Jo  (II.  10) 
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(II.  11) 


which  clearly  show  the  independence  of  t.  By  interchanging  the  order  of 
integration,  and  using  the  relation: 


T /,  \|  2 2.  1 ikxcos©  ,, 
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(II.  12) 


we  have 
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(II.  13) 


where  K - g/U^  and  in  the  last  step  terms  of  0(>^)  have  been  neglected. 
Eq.  (II.  11)  then  becomes 


poo  pit  poo  2 

✓ 1 , . _ , , , , „ -hk  + ikx  cos©  ,l7  4 ^ (sec  6)  kdk 

X.  - 5 — Urn  r(r| ) drj  d©  e cos(Ky-r))sin©>  * V- 

2r  U u-»0  ‘ k-(K,se c©  + 2ife 

r J-co  Jo  i.o  U 


-58- 


If  wc  let  j l— *0  before  we  evaluate  these  integrals,  we  then  obtain  the  same  ex- 
pression for  4 as  given  by  Eo,  (l«>)  which  is  derived  by  neglecting  the  viscosity 
from  the  beginning.  Cn  the  other  hand,  it  is  interesting  to  see  what  the  value 
of  will  be  as  x— » + cn  if  we  let  p — *0  after  the  evaluation  of  the  integral 
with  respect  to  k . First  we  rewrite  (II.  14)  as  the  real  part  of  the  following 
integral: 
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(II.  15) 


This  expression  of  4 indicates  very  clearly  that  the  surface  elevation  consists 
of  plane  wave  elements,  each  component  moving  in  a line  making  an  angle  0 
with  Ox.  Now  we  rotate  the  axes  0 x,  0 y to  Ox',  0y'  given  by 

x'  = x cos  0 + (y-'l)sinO  , y*  = -xsinO+(y-'n)cos0  (11.  16) 

as  shown  in  the  following  figure. 


In  order  to  study  the  behavior  of  4 as  x'-»  * oo,  the  integral  with  respect  to 

k may  be  transformed  by  contour  integration  considering  k to  be  complex. 

For  x'  < 0,  wc  can  deform  the  contour  to  the  negative  imaginary  axis  of  k ; 

and  for  x'  > 0,  to  the  positive  imaginary  axis.  Note  that  the  integrand  has 

■»  pu 

a simple  pole  at  k = r\  secw©  + i-y-  sec  0 which  is  located  inside  the  contour 
for  x' > 0 when  p f 0.  The  residue  of  the  integrand  at  this  pole  is  equal  to 
2 L 

-hksec  0 + ix'^.scc'O  2,.  , . 

c f',  sec  0 as  p-+0  and  should  only  concern  the  integral 

x'  0.  By  making  p zero  after  the  deformation  has  been  carried  out,  we 
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then  obtain  different  represent ations  for  x*  < 0 and  x’>  0 corresponding  to 
different  behaviors  'o*tr<  v.d  downstream  of  each  wave  element: 
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(II.  1?) 


The  first  term  in  (II.  IT)  lepresents  simple  waves  on  the  downstream  side  of 
the  wave  front  x*  - xco^P  + (y-r])  smO  = 0 ( - -j  < 0 / y ) while  the  other 

terms  represent  dist  irbanccs  symmetrical  with  respect  to  x1  = 0 and 
diminishing  exponentially  with  increasing  distance  from  x'  = 0.  It  should 
also  be  pointed  o it  litre  that  the  integration  with  respect  to  t)  (or  0)  should 
be  divided  in.o  two  regions  in  which  the  values  of  T)  correspond  to  x' > 0 
and  x' < 0 respectively  for  an  assigned  point  (x,  y).  The  final  expression 
of  K is  then 
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(II.  18) 
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whe  re 

x'  = x cos  8 - {y-r))sin0  . (II.  19) 

In  this  way  we  may  study  the  behavior  of  4 as  x-*  + go,  1 he  second  and  fourth 

integrals  with  respect  to  rj  are  negligibly  small  for  large  positive  x and 

tinite  fixed  y as  can  be  seen  from  the  fact  that  both  upper  and  lower  limits  of 

integration  tend  to  + oo  and  -oo  respectively.  The  same  is  true  for  terms 
. , -x'q  -ic' q 

with  c and  e which  may  also  be  neglected.  Hence  the  only  significant 

term  in  Eq.  (II.  18)  for  large  positive  x is 
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(II.  20) 


which  is  the  same  as  Eq.  (20).  The  limiting  value  of  the  second  integral  goes 

to  zero  as  0(  -rr)  which  c in  be  shown  by  applying  t lie  method  )l  stationary  phase, 
/x 


* = * 


oo 


sin( 
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rln)  an 
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as  x ->  + oo  , 

and  y fixed,  finite 


(II.  21) 

For  x negative  and  large  with  y again  finite  and  fixed,  the  first  and 
third  integrals  are  negligibly  small.  The  sum  of  the  other  two  integrals  can 
be  shown  by  using  Watson's  lemma  to  be  much  less  than  the  absolute  value  of 
4 given  by  Eq.  (II.  21). 


It  is  of  interest  to  note  that  by  considering  a convergence  factor  ^ 

(usually  called  Rayleigh's  viscous  term)  in  the  intermediate  stages,  the  final 

result  of  surface  elevation  becomes  asymmetrical  with  respect  to  x : 

4 damps  out  like  0 (- y)  on  the  upstream  side  and  has  a wave  formation  down 

stream  of  the  hydrofoil  and  finally  diminishes  like  0 (-4*)  for  x large. 

jx 
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APPENDIX  III 

Calculation  of  Wave  Resistance  by 
Method  of  Traveling  Pressure 


According  to  the  method  of  traveling  surface  pressure  (cf.  Roy.  Soc. 
Proc.  A,  Vol.  93,  p.  2-14  (1917)  ),  the  wave  resistance  is  simply  the  total 
resolved  surface  pressure  in  the  x-direction.  Because  4 is  everywhere 
small,  this  leads  to 

R = H F(ffl)  ^dx  dy  (III.  1) 

taken  over  the  whole  surface.  F ( cd)  is  given  by  (II.  2)  to  be 
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(III.  2) 


f 2 2 ’ 1 3/2 

Since  j x + (y-h)  + h~  | is  even  in  x and  (y-q),  »t  follows  that  only  the 

part  of  3$.  even  in  both  x and  (y-T})  contributes  to  R.  Using  the  form 
dx 

of  4 given  by  Eq.  (18),  which  is  the  only  part  even  in  x and  (y-q),  we  have 
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Now  the  integrals  with  respect  to  x and  y can  be  carried  out  by  introducing 
the  transformation: 

x = r cos  0 y - rj  = r sin0 

and  using  the  notation 

a = r(  sec  0 b = /.  scc^O  sinfl 
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wnere  f and  g are  defined  by  Eq.  (27).  This  result  agrees  with  Eq.  (29). 
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APPENDIX  IV 

Evaluation  of  Some  Integrals 

(A)  Derivation  of  the  Integral  in  Eq.  (46) 

Substituting  Eq.  (45)  in  (43),  integrating  by  parts  with  respect  to  t} 
and  using  the  relations 

y = b cos  0 rj  = b cos  f , (IV.  1) 


we  have 
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The  interchanging  of  the  order  of  51  and 
resulting  scries  converges.  Noting  that 


sign  is  justified  if  the 


r w 


sin(z  cos  ) cos(2n+  1)  jdi)r  = (-)n*  ^2n+  1^Z^ 


(IV.  2) 


where  ^n+l  ^cnolos  Bessel  function  of  the  first  kind  (cf.  Ref.  24,  p.  20), 
we  finally  obtain 
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which  is  Eq.  (46). 
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(B)  Derivation  of  the  Integral  in  Eq.  (4?) 

Substituting  Eq.  (45)  into  (44)  and  proceeding  in  a way  similar  to  that  in 
the  previous  case,  \*»c  have 
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The  first  integral  on  the  right  hand  side  can  be  integrated  to  give  (cf.  Ref.  25, 
p.  37  and  p.  33): 
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Combining  Eqs.  (IV.  4-6),  we  obtain  Eq.  (47). 
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(C)  Evaluation  of  the  Integral  in  Eq.  (61) 
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where  K and  c.  denote  the  complete  elliptic  integral  of  the  first  and  second 
kind  respectively. 


(D)  Derivation  of  the  Result  of  Eq.  (64) 


Applying  the  transformation  sec  9 = cosh  ^ to  the  integral  in  Eq.  (63), 
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Note  that  the  above  integral  still  converges  as  — »0  , keeping  6 constant, 

^ •) 
One  way  to  obtain  the  result  given  in  Eq.  (64)  is  by  expanding  the  term  Jp 

into  an  infinite  series  (cf.  Ref,  24,  p.  147)  as  follows: 
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lntroducing  this  expansion  to  Kq.  (IV.  8)  and  interchanging  the  order  of 
an.l  Tj  sign,  we  obtain 
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provided  that  the  resulting  series  converges  absolutely.  The  region  of  conver- 
gence for  a and  \ will  be  shown  below.  By  using  the  following  relations 
(of.  Ref.  24.  p.  172) 
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where 


K denotes  modified  Bessel  function  of  the  second  kind,  we  obtain 
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It  can  be  shown  that  the  above  series  has  different  regions  of  convergence  for 
large  and  small  values  of  a due  to  the  different  behavior  of  the  function 
Kn(~).  When  <r  is  small,  this  series  converges  absolutely  and  uniformly  for 
any  close  interval  of  (3  ; but  for  large  values  of  a ( >>1),  the  above  series 
only  converges  for  £ - k^I-  The  interchanging  of  the  j and  ^ opera- 
tion is  justified  only  when  the  above  conditions  are  fulfilled. 


(L)  Derivation  of  the  Result  of  Eq.  (68) 

We  want  to  sum  the  series  in  Eq,  (IV,  10)  for  large  values  of  <r.  In 
this  case  K (— ) has  the  following  expansion: 
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Ko(i)  = log  (2<r ) yi)*  £ i&lf™  v ( m+  1 ) , 

m=o  (mW 


y 


(IV.  11) 


+ ^(m+l)  + 2 >jr  (n+m-fl)  | , (n  > 1) 

It  is  easy  to  see  that  the  most  important  contribution  to  the  sum  of  the  series 

in  Eq.  (IV.  10)  comes  from  the  term  y(nl  )( Zvf1*  * of  the  expansion  of 

K .,(•*•).  The  detail  of  the  calculation  ran  be  shown  as  follows: 
nr  i a 


Decompose  V into  three  parts  such  that 

J<  - 4 \ 4 

i -^11  -^12  • 


o 


where  and  -‘Jj  3 ar0  given  below. 


..  e"?  " »-fV<n+T)|  W>  , " f 1 , 1 


(IV.  12) 


e «■  ^ (“Ff  (n+4);  (n+‘|’)  ^ n+1 

= 47  ^;nl"(n+l)(‘:(n+2)l  ^ 

n=o,l 


! nl  (2<rf+1  i 


c » 1 (-fri( n+,|)  j n+ 1 | (-Pp2(n+-|)  1 n 

TT  } /-  In+H'In^)!  ^ * “TZ_  n(nVr)rrn^)l  ^ 

*•  n-o  * ’P  n=l  * 


l6p 


log(2<r ) i 


(IV.  13) 


Now  both  these  two  series  converge  only  when  X.  y 1 ; but  the  first  series  can 
be  directly  related  to  a liypergeometric  function  which  may  then  be  continued 
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analytically  to  the  region 

<&-  (-f 1 2(n+4  i n+i 

2-  Tn"TTlF+"2): 

n=o  K 


- TT 


1 - 


0 < X 1 as  follows: 


2/ 1 


(?)  ^M-fr  *<•♦*) 

r(2)  f-  r (.+2) 


a =o 


v/1  1 • ? • 1 i 1 

A (7'  7 * 2 * *77' 


x rfl  3 - - • 1 

1 r VT  . -T  • 2 . 

2 4 2 1+X' 

•l+X  1 * 


1 


If  we  use  one  of  Gauss's  recursion  formulas,  the  above  hypergeometric  function 
can  be  expressed  in  terms  of  complete  elliptic  integrals,  (cf.  Ref.  26,  pp.  9,  10). 


,1  3 


,.,1  1 


1 1 


,2.  1 


^ (7 . 7 ’ 2 • k ) = “I  * ^7  * 7 ; 1 ; k ) j 


= 1 k'2  K(k)  - z(k)] 


F inally  we  have 


00  (-] P I "2(n+-j)  j n+l 

>■  (n+ 1 JT  (n+7)'.  (“7> 


n=o 


) ’• 


- x i 1+x2 

IT 


K ( ;=—)  - E ( 

i 1+X2  Jl+X2 


(IV.  14) 


To  carry  out  the  summation  of  the  second  series  in  Eq.  (IV.  13),  we  note  that 
if  wj  write 


f(z) 


n(n+lj!  (n+2)'. 

n=l 


-n 

z 


z - X 


(IV.  15) 


then  f'(z)  may  be  expressed  in  terms  of  a hypergeometric  function  together 
with  an  algebraic  function  as  follows: 


(■(.)  = 


co  (.frV4)  r 

2.  (n+ 1 )! Tn+TTT  z = ' T 


n = l 


2/ 1 » _ 

'7*  | . 11  r,l  l.  2.  U 
1 - I • r'7*  7’  z'  V 


1“  (2)  L 


7 


- tr  1 - J-  - i r*(i+*)  [ K(-iT-)  - E(-4r-)  i i . 
I 8 Jl+z  A+z  J j 


(IV.  16) 
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As  z = X2 -*  0,  . -*  1 ; the  asymptotic  value  of  f'(z)  as  z-*-0  may  be 

, 1 + z 

deduced  from  the  above  relation  by  using  the  known  asymptotic  expansions  of 
K and  E. 


*'(*)  = IT  - 1 + 


SI 


+ - 

IT 


z_ 

+ z 


log  4 


Tiz 

z 


+ 0(z  log  z) 


Since  the  integral  of  an  asymptotic  expression  is  also  asymptotic,  the  asymp 
totic  value  of  <(z)  can  then  be  obtained  by  integrating  the  above  equation. 

The  final  result  is 


f(z)  = + £ log  z + 0(  1 ) = + £ log  £ + 0(  1 ) 


(IV.  17) 


Therefore,  substituting  Eqs.  (IV.  14,  17)  into  (IV.  13),  we  obtain: 

^11=  \ i 1 |K(— =r)  - £(^=r)|  * -Ulogp+O^,  i)) 

1 fi+x2  W 1 r aj 

The  second  part  of  , namely,  ^\i  • *s  8*ven  Py 


t — e 
e>  u-- 


a oo  (.ft  (n+  j)Rn+|)  ^ ^ n 

TTI  2_,  nl  (n+  i j:  (nt  2)1 
n = o 


k„+i<7> 


(nl  ) (2<r) 


n+  1 'i 


eT  ® (-fr(n+-|)I  (n+|)  a n+ 1 

717  2-  " nl  (n+  l )l  (n+  2)T ""  {TT) 
n=o  <-P 

e*7  °°  (-fr(»-«j)r(*+-j)  ('7) 

i!2) 


4»r<r 


e 

’ 77 

__i_ 

<r 

e 

•7T 


V 

s = o 


(s+2) 


1 TT  / ^ 1.1.  1 \ 1 

l-F(-7,  7.  2.  --j)  j 


jl  + X2F(-i,  -r;  2;  — j ) - X 

L C 1+X4 


(IV.  19) 
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After  this  analytic  continuation,  if  one  use9  again  Gauss's  recursion 
formula,  one  finds, 

1 3 , 1 . _ 2 .2  I , 1 3 . , . 1_.  1 3 1 x \ 

f(-4i  7 i 2 ; j)  - t\  ^ F(t,  T • 1 • — ~T/  ‘ F ( * 7 ■ 7 • 1 • ' 

1 1 1+XT  3 l L L l+\  1+X  j 

= 4 £(itn2)  e( — — ) - X2  F ( - i . 4;  1 i-i-r)!  (IV.  20) 

T i " rr<7  i 

Combining  Eqs.  (IV.  19)  and  (IV.  20),  we  have 

J,2  * \ - 4x  - x2  J nx2  rt-f  fl  |(l+o(i)) 

V ' 

(IV.  21) 

It  is  evident  that  the  hypergeon  etnc  function  in  the  above  equation  has  singu- 
larity at  X.  - 0.  To  calculate  its  value  near  X - 0,  we  need  the  following 


formula 


oo 

r'(a,  b;  a+b;  z)  = an(l-zf  P„  + log(l-z) 


(IV.  22) 


whe  r c 


1 (atb)i  (a+n)|  (b»n) 
P (a)  J(b)r2(l  + n) 


Pn  = (a+n)  + •;  (b+n)  - 2 V (1  + n)  . 


Hence 

1 3 , 1 . 2 j , , . 4,'l  + X^.  3 ,1  , 4 /T+X2  w X^_v  1 . n/v4  . 

F(-y,  7 ; 1 ; t)  = - (2  - log  — L ) - j-  (t  - »°g  — ^ )( O +0(X  logX), 

L L l+X^  51  X 3 X 1 + X* 


and  finally  we  have 


(1  n 

The  third  part  of  denoted  by  Oxj3»  *s 


(IV.  23) 


Jl2=^  ((nx^-^x-’x2  W «og  *]$£  ♦ o<x2)  f-  (iv.  24) 
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( _ 
v S 


7 ® (-/TV -1)1  (n+4) 


n 


13  2 -T^X)'' 

4tip  n = o 

e‘ff  [ (-f  f(n4  i-)[~(n+-|) 


(~~T ) K (±) 


20- 


no" 


V II 


(-V)"  ,(n-l)l(2j) 


+ £ log  2<r 


or, 

_ _1_ 

- oo 

f 

e I 

1 n^l 

o*  13 

*Ha. 

t-. 

CO 

The 

series  in  the 

that  applied  to  Eq.  (IV.  15).  The  final  result  is 


/; 

^13 


16 


1 logp+  0(-l  i)  . 


7 


(IV.  25) 


P' 


Therefore,  substituting  Eqs.  (IV.  18),  (IV.  21)  and  (IV.  25)  into  Eq.  (IV.  12), 
we  obtain  a good  approximation  of  the  value  of  0>  ^ defined  by  (IV.  8)  or 
(IV.  10)  for  a large  as  follows: 


•^1™ 


~jt/  Z 2 2_ 

< — sec  © 

<r 
e 


; s(Tg  sin  8 ] 

-L.  ' 

.••»«.  c 0 sin  6 


sec  © d © 


-J.  ( . - Ju\-  ; K(  ) . i(-^r)  1 [ +^||<HX2) 


- iJl-i 
l * 


l+x 


J l+\2  ' J 


J-2 


1 + X. 


_ 4 X - X2  ; U\2  F(-i,  y ; 1 ; -Lr)  ' + — U logP  + 0(iy.  ±)  (IV.26) 

Z Z l+X^  64  p^  pr  ff 


(F)  Derivation  of  the  Result  of  Eq.  (73) 

Using  the  result  given  in  Eq.  (IV.  7),  we  may  write  the  integral  in 
Eq.  (71)  as  follows: 


-72- 


AL, 


p'o 


2 f tr/2  | .ff/2 

csc28 a#  ■!  l - i — £2i 


l "l 


sin2 8 sin*0  J 


p r; 


it/ 2 


cot  6 


1-i 

TT 


COS 


vn/2  rit/2 

. 4 


j 1 + -^  sin2®  sin0; 


ttX 


cos  0 dO 


8-0 


ain20  cos20  d0 


iQ  | l+-i>  sin^8  sin20^ 


3/2 


2 fir/ 2 

^ ! o 2n  Jn 

— -j—  cos  8 d8 

tt\ 


r/1 


sin2t(  coa20  dfl 


.1  2.  1 . 2Q  . 2d  > 

it-  sin  fl  - A sin  Osin  ^ 

»** 

\ X. 


3/2 


PI 


ttX 


r/2 


cos2  8 


6 { — de 


sm28  J3/2  ^f7Tsin20 


(IV.  26) 


where  <£(k)  ’S  a derived  complete  elliptic  integral  (cf.  Ref.  25,  p,  73) 
defined  by 


6(k)  = 


tt/2 


La  La 

sin  0 cos  0 d0 

f.  .2  . ZA  3/2 
| I - k sin  0 i 


(IV.  27) 


finally,  changing  the  variable  of  integration  from  6 to  k by 


k - 


sin  8 


. X2  + sin2 8 


(IV.  27a) 


we  obtain 
AL 


1 


I 1 - £)2 

| T-  <*(k)  dk 

1 - k‘ 


(IV.  28) 


whe  re 


r = 
y* 


ii+x2 


(IV.  28a) 


This  is  the  form  used  in  Eq.  (73), 
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(G)  Derivation  of  the  Result  of  Eq.  (76) 

Substituting  the  expansion  of  6(k)  given  by  Eq.  (75)  into  Eq.  (IV.  28), 
then  applying  the  transformation  k2  = X2t  and  integrating  termwise,  we  have 


t + 


75  , £2  2.  ... 

T (T}  1 + 245 


= - »Prc 

o4k 


o T 


2F1 


4' 


1 

7' 


2. 

) 


( 2 

4> 


2F1 


4 


3 

V 


f- ) 


..2  4 


(?•  7;5  :^2)+0  ((£)  j 


I 


(IV.  29) 


where  (a,  b,  c;z)  denotes  the  hypergeometric  function  which  has  the 

following  integral  representation: 


p r 1 

^(a.b.c;*)  = |'(b)^.b)j  tb~  1(l-t)C"b"1(l-2it)_adt.  (Rec>R2b>0). 


Expanding  these  hype rgeorae trie  functions  into  convergent  series,  we 
finally  have 

_ n(>ro  i , . 5 1 181  1 3015  1 0f  1 }\ 

wJr 

(IV.  30) 


(H)  Evaluation  of  the  Integral  in  Eq.  (79) 

Applying  the  transformation  secO  = cosh  ^ to  the  inner  integral  in 
Eq,  (79),  and  proceeding  in  a similar  way  as  shown  in  (IV.  D),  we  have 
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/•n/2  2u  2_ 

! - — sec  0 


^2 


(u)  = 


,u  2 


sec  0 sin  0) 

2 

•v  sec  V sin  0 


sec  0 dO 


1 <r  j 

7 e 


U r OO  U . 

1 - — cosh  u 


\ sinh  u) ) 


sinh  u 


, 1 + cosh  u.  . u , 

( 2 ~ - ) cosh-j  QU 


I (-rnn^)  (Ui2n*2 

e n'.(n+l)!(nTZp.  ^ 


OO  u . 

- — coshu  y 

e 9 sinh  nu(l+cosh  u)  cosh^  du 


n=o 


Now  sinh^nu  can  be  expanded  in  terms  of  cosh  mu  as  follows 


2n 

. ,2n  (e  -e  ) 
sinh  u = » — y * '"  = 

2 


n- 1 


y—  — = -y-  ; 2n  ) ♦ — >,  V (-fn  ( 2n  ) cosh(2n-2m)  u 

**  \ n ! \m  / 


m=o 


and 


( 1 + cosh  u)  cosh  ^ = -j  cosh  ^ cosh  . 

Using  the  relation  (cf.  Ref.  24,  p.  181) 


foo 

-z  coshu  , , 

i e cosh  v u d u = 


Kv(’i. 


(R i .>0,  Rf  i>0)  , (IV.  31) 


we  then  have 


n- 1 


. U 00 

( 

O’  N 

THT 

n=o 

(-r 

/ 2n  ' 

V rrV 

* * A V 

1 1 •’ 

* E <-*"  vt) ! 3K2(„-"<)ti<7> ♦ 3 ‘V-i- E?' 


m=o 


+ K2(n-m)-i(;>;: 


2 

(IV.  32) 


When  c is  very  large,  then  (cf.  Ref.  24,  p,  80) 


u 


Knti<T>=  ]%•'  <>♦“#  /z?'’7  <lv-»> 

2 '7* 
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It  can  be  seen  that  the  first  order  term  of  ^^(u)  comes  from  the  term  with 
K^n+  3 (j).  hence  for  large  ff,  J? ^(u)  may  be  approximated  as 
2 


l~n~  "Tl  (-fr(n+4)R2n+v)  ,^2n+2 
, irr  e i L ~nrrKVTm^T2'H~'  w 


2^u^  = IT^  e 


2u 

i I I u <r 

ff 


T«  J ~ c 


^ n=o 

( ^>) 


nTTn+TTHn+TJ!  (^} 


0+  <H±)) 


1-  p^:  F^.l;2;-^) 

l^2  l + 5?. 


O+0(£)  ) 


. Ri)  nr  -¥  r , r<|) 

4 fir  J * l rfRl 


. ,==!  (I  + 0(X2  i)),  ff  -*  oo 

Ji+v2i 


(IV.  34) 


For  small  values  of  a (h,  X,  large),  it  can  be  shown  by  the  method  of 
steepest  descent  that  has  the  same  asymptotic  value  as  ji  ^ given  in 

Eq.  (66),  that  is 


2u 


3 

7 


J»S  ‘ e'  ' /if  2.  (,. 


z'"1  ~ ^ w l ~~T  — 2 ( i + o(ff)  ) * ttS  • 


(IV.  35) 


Therefore,  for  a large,  we  have  from  Eqs.  (79)  and  (IV.  34) 

2u 


4 j it  [7  n f (^)  jt+x  ■. 


r<T) 


00 


— du  ( 1 « 0(X2,  i)  ) 

fu(u-l) 


= id-  -r^r  ^ 

4 Jir  j«r  Jwj  (7)  ;i+X  ' 


r<*> 


2 00 


I n<r  a 

* J ~T  +e 


Sle^'dt 


) 

(1  + 0(X2,  i)  ) 


pfor  i 1 r(|) 


4 J w 4 i«r  J «["(?)  /l+x2  l 


2 

fwff  . <r  ^ #2v 

J T * ' £-.(7) 


(ito(x2,±)). 


where  E.(t)  represents  the  exponential  integral  function  defined  by 


£t(0  = 


00 


-X 

e dx 

x 


(t  real,  cf.  Ref.  23,  p.  471) 


-t 
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w'nich  lias  the  following  expansion 


;,<7>  * 


+ log 


2 . 2 . 1 ,2. 2 . 1 ,1?  . 

7 + 7 + ITT r (7>  + 7HT  (7}  + * * 


, , for  » large. 

(IV.  36) 


Finally  wc  have 

Pf"o  - 1 I'(t) 

*4  =— s-i  4)(>-— ==ri 

4J2  J"  <?>  ,W; 


K fe<r  + i°8|> 


( 1 ♦ 0(  \2,  i)  ),  as  <r-*oo 
(IV.  37) 


*4 


On  the  other  hand,  for  <r  small,  the  asymptotic  value  of  aL^  is 

_ 2 . uo  2 u 
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O TT(T  I j tT  y u 


pr„ 


JT  ~T 

P 


^—j-du  ( 1 + 0(ff ) ) 

1/2 


rsr  i ‘I  f * <*  + f*>  dt  (lt0(  .. 

— e | e j dt  (l  + 0(r)) 

p J.l 

a 


pr0 


j e 1 “i<;>  ( i + °(<r) ) . 


The  asymptotic  expansion  of  Z.(^)  lor  J-  large  is 


2 

<r  /<r 


E.(f)  = e"  (J)  (l+f+  21  (|)  ♦ 0(0) 


Hence, 


r 


,2 


*4  = 


3 

7 5 , 

J v (^)  (i  + 7 + ) )* 


as  <r  0. 


(IV.  38) 


(I)  Derivation  of  the  Result  of  £q.  (83) 
From  Eq.  (82a),  we  have 
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4\ 


Now  if  we  continue  the  hype rgeome trie  function  analytically  to  region 
OCX  Coo,  we  obtain 


Uj(y)  = - 


r. 


i 


4 fir  b , . 1V2 

2X  J 1 + 4X  m =o 


v—  (*rh)‘'nRm+|)  j m i i i 

L m*  1 + 4X^  L C 1+4X 


1+4X 
(IV.  39) 


(J)  Derivation  of  the  Result  of  Eq.  (85) 

Applying  the  transformation  ^ = (bCsec^O)u  to  Eq.  (82b),  we  have 


/ \ ft  ' o 

UA> 2 • — 


where 


,oo 


du 
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(IV.  40) 


O’  = 


TT7T  ' 


p = rz  • n=  f • 
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The  inner  integral  can  be  treated  in  a way  similar  to  that  applied  to  in 

(IV,  H).  In  particular,  at  y = 0,  we  have 
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For  <r  very  large,  the  first  order  term  of  u2(0)  comes  from 
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oince 


<T*r<7> 


‘1  sin 


TrT 


J2  IT 


F(l’T;2;--iz)  ~—^x)r!/i  F<f  4: 2;-i-r> 

4X  ( l+^V  )^7  4 1 1 + 4X4 

and  (cf.  Eq.  IV.  37) 
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(Y+  l°g  “) 
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we  finally  obtain 


u2(0)^  ‘ TE  TT/l  F^T'  T’2’ — 

2 4b  J 2X  (1  + 4X2)V  4 4 1+4X2 


3 3 


M(y+>o8t)  ! <i+°(?)) 


J 


(IV. 41) 


This  formula  is  good  for  all  X,  O^X^oo  , provided  that  <x  is  large.  The 
above  hype rgcometric  function  can  be  expressed  in  terms  of  tabulated  elliptic 
integrals  by  using  Kummcr's  transformation  and  Gauss'  recursion  formula  as 
follows: 


r ,3  . 3 . _ , 1 y _ ,3  3,1 

f(T’T'2, V ' * '7*  t • 1 *t  * ; 

1+4  x 


1 +4X 


(Kummcr) 


2 1 r/1  3 . . 2 . _,1  1 . .2. 

= ~7  F(7’  V 1 ;k  ) ' F(7  * 7:  1 ;k  > 

k 


(Gauss) 


whe  re 


k = 


1 X 

~ 

j 1 + 4X 


Applying  the  transformation  formula  once  more,  we  obtain 


w3  3 , 1 N 2 f 1 „,1  1 , ,2X  ,.,1  1 , ,2, 

T'  , T?  7’ ' ; 7’  ' 


4 1 

TT  77 


k (1_L 

4 2 J1+4X2 

“ f - 

J 1+4X  + 2X 


E(k)-(l-k2)K(k)l  =| 

1 -K 


B(k) 


B 


» 1 

7 


/ ■> 

V 1 + 4X 
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whcre  B(k)  is  defined  by  £q.  (7  in).  Substituting  this  relation  into  (IV.  40), 
we  hive 


2 ; 


o 1 


l 


l 


^ ITx  -+4xi  + 


- B ( / i - -===.  N(l+0(i)  ) . 

I I 2 I 2 1 a 

Z\  W j 1+4X*  ' 

(IV,  42) 


When  <r  is  very  small,  we  can  apply  the  method  of  steepest  descent  to 
the  inner  integral  of  £q,  (IV.  40),  The  result  is  as  follows: 

ao  2u 

« a du  (1  + 0(<r)  ) 


(cf.  Section  H,  £q.  IV.  38) 

(IV. 43) 
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J Zn<r 
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(K)  Derivation  of  the  Result  of  Eq.  (91) 
The  method  used  to  approximate 


jt/2  2 Zn  . ,1  2_ 

i sec  0 J,(»sec  OsinO)  c 

♦ ^(0,0,-h)  --  -;.r  ; -S2, 

o 


sec  P dO 


sec  8 sin  8 


lor  <r  large  is  very  similar  to  that  discussed  in  (IV.  D and  E).  Introducing 
the  transformation  sec  8 = cosli  ^ to  the  above  integral  and  then  expand  Jj 
into  a series,  integrating  termwise,  we  obtain 


9^(0,  0,  -h) 


r - — l03  - -coshu  J,(ttt  sinhu)  . . .2 
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8 frr  p 


i 00 
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r — . ( - ! ( n+  -r)  n * . i i 

^ ’ ntCnTTT T {~h)  Kn(7)+  ( ?) ) Kn+ 1(7} 

n-o  op  >. 


(IV.  44) 
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where  the  relations  given  by  Eq,  (IV,  9)  have  been  used.  Referring  to 
the  expression  of  K (— ) for  <r  large  given  by  Eq,  (IV,  11),  it  is  easy  to  see 

**  I j 

that  the  first  order  term  of  comes  from  the  term  y(nl  )(2<r)  of 

1 L 

the  expansion  of  K .,(—).  Then 
r n+l'<r 


(0.0,h)=  - Ilia 
* » n<ft 


■7  “f  , ,1.»U  P 

c ^ "nrtn+rjr  {7Z]  r(l n+7>7(2<r> 


i+i  ii 


i<r0  . Vs  i "+1 


2-  fn  n = 


« HUn+T,  n+l  r , 

P ^ (nil)!  Z>  l 1 + > J 


1 gA-V}  i> 

i nTo  r (4)  nl  J I P 


i - •!}  - i ( i * o(i,  1) 

' WJ  ' ' ? 


(IV. 45) 


where  the  appropriate  analytic  continuation  from  to  O'lX^oo  is  made 

at  the  last  step. 
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